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THE  VERIFICATION  OF  A  COMPUTER  MODEL  OF  INTERNAL  LIGHT  REFLECTIONS 


FOR  HELICOPTER  CANOPY  DESIGN 


INTRODUCTION 

The  US  Army  Human  Engineering  Laboratory  (USAHEL)  has  developed 
several  computer  programs  for  evaluating  helicopter  canopy  designs.  The 
programs  provide  a  point-wise  computation  of  the  internal  glare 
reflections,  external  solar  glint  and  the  optical  effectiveness  of  the 
transparent  surfaces  of  a  helicopter  canopy  design.  The  programs  were 
developed  using  geometrical  ray  tracing  methods. 

The  purpose  of  this  report  is  to  verify  the  application  of  ray  tracing 
concepts  to  the  computation  of  canopy  performance.  The  report  describes 
the  method  used  to  test  a  computer  model  for  the  internal  light 
reflections  on  the  transparent  surfaces  of  a  helicopter  canopy  and  the 
test  results. 


METHOD 

A  mockup  of  the  Model  209  AH-1S  Helicopter,  flat  plate  canopy  design, 
was  placed  in  the  enclosed  bay  of  a  building.  A  coordinate  reference 
system  was  laid  out  on  the  floor  about  the  moocup.  An  observer  sitting  in 
the  pilot's  position  of  the  mockup  used  a  pointer  to  pinpoint  the 
locations  of  direct  view  and  reflected  images,  on  the  canopy  surfaces,  of 
an  outside  light  source.  In  the  first  phase  of  this  test,  the  luminances 
of  the  images  were  measured  with  a  telephotometer.  In  the  second  phase, 
the  coordinates  of  the  locations  of  the  images  on  the  canopy  surface  were 
measured.  The  data  were  measured  at  different  positions  of  the  light 
source  on  the  outside  floor.  The  source  positions  were  different  for  the 
two  phases. 

The  transmittance  and  coordinates  of  the  image  locations  were  computed 
for  the  different  light  source  positions  using  a  computer  program.  The 
measured  and  computed  values  were  compared  by  a  student’s  t-test  for 
matched  samples  and  by  a  linear  regression  analysis. 


TEST  FACILITIES 

The  helicopter  mockup  was  supplied  to  USAHEL  by  the  Project  Manager's 
Office,  USA  Aircraft  Survivability  Equipment  (ASE).  The  mockup  had  been 
built  in  accordance  with  the  specifications  for  the  Model  209  AH-1S 
Helicopter,  flat  plate  canopy  design.  The  mockup  was  placed  in  an  enclosed 
bay  area  for  the  test.  The  room  was  sealed  for  light  control  and  the  walls 


and  ceiling  were  painted  with  a  flatbaae  black  paint  to  prevent  extraneous 
reflections.  The  two  ceiling  lights  were  shielded  to  prevent  direct  glare 
and  the  lights  were  dimmed  during  the  test. 

A  brow  pad  was  Installed  in  the  cockpit  of  the  mockup  at  the  pilot's 
position.  The  pad  extended  vertically  downward  from  the  top  canopy  surface 
and  was  positioned  so  that  the  eyes  of  an  observer,  sitting  in  the  pilot's 
seat,  would  be  in  the  nominal  position  when  his  brow  was  against  the  pad. 

A  rectangular  (x,y)  coordinate  system  was  laid  out  on  the  bay  area 
floor  for  positioning  the  light  source  during  the  test.  The  origin  of  the 
coordinate  system  was  positioned  directly  under  the  nose  of  the  mockup 
with  a  plumb  bob.  The  y  axis  was  laid  out  directly  beneath  the  centerline 
of  the  mockup  parallel  to  its  longitudinal  axis.  The  x  axis  was  laid  out 
normal  to  the  y  axis.  The  x  and  y  axes  were  marked  in  1-foot  intervals  for 
use  when  positioning  the  light  source  during  the  test. 

A  5-foot  stand  was  built  with  a  movable  clamp  to  hold  the  1.5-inch 
diameter  light  bulb.  The  stand  could  be  moved  to  any  position  on  the  bay 
floor.  The  clamp  holding  the  light  bulb  to  the  stand  could  be  raised  or 
lowered  to  any  height  up  to  5  feet  above  the  floor  level.  The  stand  and 
the  light  bulb  served  as  the  light  source  for  the  test. 


Apparatus 

The  apparatus  employed  in  this  test  were  used  to  measure  (1)  light 
luminance,  (2)  the  coordinates  of  the  canopy  surface  points,  and  (3)  the 
coordinates  of  the  vertices  of  the  mockup  frame.  The  apparatus  and  their 
usage  are  as  follows: 

1.  Light  Luminance 

a.  The  luminance  of  the  light  source  images  on  the  canopy 
surfaces  was  measured  with  a  telephotometer  with  the  following  components: 
Model  220-31  photometric  telescope  (aperture  setting,  20  minutes),  Model 
820A  digital  photometer.  Model  820-18  control  module,  Model  820-18 
photoplc  correction  filter.  Model  820-18-2(S4)  photosurface  tube,  Model 
820-19  general-purpose  housing,  and  Model  700-3B  fiber  optics  probe  (Gamma 
Scientific  Co.,  Burbank,  CA).  The  instrument  error  is  5  percent. 

b.  The  telescope  was  focused  and  centered  on  the  image  of  the 
light  source  by  the  observer  who  was  sitting  in  the  pilot's  seat.  The 
luminance  was  read  from  the  digital  meter  and  recorded  by  an  assistant 
while  the  observer  kept  the  telescope  centered  on  the  image. 

c.  The  luminance  of  the  1-1/2-inch  diameter  light  source  was 
measured  from  outside  the  mockup.  This  direct  line  of  sight  measurement 
served  as  a  baseline  for  transmittance  computations  and  was  rechecked 
during  the  test. 
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2.  Coordinates  of  Canopy  Surface  Points 


The  three-dimensional  Cartesian  coordinates  of  the  light  source 
on  the  canopy  surfaces  were  measured  using  a  plumb  bob,  a  level  and  a  tape 
measure.  The  observer  sitting  in  the  pilot's  seat  used  an  expandable  metal 
pole  during  the  test  to  pinpoint  the  location  of  the  light  image  on  the 
canopy  surface.  The  position  was  marked  with  a  small  tab  of  adhesive  tape. 
The  distances  separating  the  image  position  from  a  convenient  vertex  on 
the  mockup  frame  were  measured  using  the  level,  the  plumb  bob,  and  the 
tape  measure.  The  coordinates  of  the  image  position  are  the  distances 
along  the  three-dimensional  axes  of  the  Cartesian  coordinate  system  for 
the  mockup,  which  separate  the  image  from  the  vertex. 

3.  Coordinates  of  the  Vertices  of  the  Mockup  Frame 

a.  The  coordinates  of  the  vertices  of  the  mockup's  internal 
structure  and  canopy  frame  were  measured  using  a  surveyor's  transit 
theodolite  (Keuffel  and  Esser  Co.,  Hoboken,  NJ).  The  internal  structure 
Included  the  copilot's  seat,  gun  sight  and  display  panel,  and  the  pilot's 
side  armor  and  control  panel.  The  canopy  frame  included  the  beam  structure 
over  the  cockpit  which  held  the  transparent  surfaces  in  place. 

b.  The  coordinates  of  these  vertices  were  measured  in  the 
three-dimensional  (x,y,z)  coordinate  system  of  the  mockup.  The  y  axis  was 
along  the  longitudinal  axis  of  the  mockup  with  the  positive  direction 
toward  the  increasing  station  lines.  The  z  axis  was  the  vertical  axis 
directed  along  increasing  water  marks.  The  x  axis  was  orthogonal  to  the 
other  two  axes. 


c.  The  coordinates  were  measured  using  the  theodolite,  a  steel 
rod  and  a  tape  measure.  The  theodolite  was  positioned  in  front  of  the 
mockup  and  two  reference  directions  were  laid  out  on  the  bay  area  floor. 
One  direction  was  along  the  longitudinal  y  axis  of  the  mockup,  and  the 
other  orthogonal  to  it  (x  axis).  The  theodolite  was  then  lined  up  with  the 
appropriate  axis  for  the  coordinate  to  be  measured;  i.e.,  along  the  y  axis 
for  the  x  coordinate  values,  along  the  x  axis  for  the  y  coordinate  values, 
and  leveled  with  the  horizon  for  the  z  coordinate  values.  The  coordinate 
values  were  measured  by  placing  the  point  of  the  rod  at  the  vertex  being 
measured.  The  rod  was  aligned  with  the  coordinate  axis  using  the 
crosslines  of  the  theodolite.  The  distance  along  the  rod  from  the  vertex 
to  where  the  crosslines  crossed  the  image  of  the  rod  was  measured  as  the 
coordinate  value. 


Observer 


The  test  is  of  a  mathematical  model  for  light  reflections  and  no 
subjects  were  employed.  The  single  observer  was  one  of  the  data  collectors 
for  the  test.  His  ability  to  place  a  pointer  at  the  position  where  he 
observed  a  light  image  on  a  surface  was  not  considered  to  be  a  test 
variable. 
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Procedure 


The  test  was  divided  into  two  phases.  In  the  first  phase,  the 
luminance  was  measured  for  all  direct  view  and  reflected  light  images.  In 
the  second  phase,  the  coordinates  of  the  canopy  surface  positions  were 
measured  for  the  direct  view  and  reflected  light  Images.  The  two  phases 
used  different  sets  of  light  source  positions. 

A  set  of  15  test  positions  was  selected  at  random  for  the  light  source 
used  in  each  test  phase.  The  set  was  Selected  from  a  larger  set  balanced 
with  respect  to  floor  location  and  stand  height  about  the  mockup.  The  set 
was  also  balanced,  as  much  as  possible,  with  respect  to  the  windows  and 
canopy  surface  locations  where  reflections  were  expected  to  occur.  The 
test  set  was  randomly  ordered  for  presentation  to  the  test  observer. 

The  light  source  was  presented  to  the  test  observer  from  the  selected 
test  positions  in  a  series  of  test  runs.  The  test  assistant  positioned  the 
light  source  at  the  floor  location  and  stand  height  specified  for  the  test 
run.  The  stand  was  located  on  the  floor  by  tape  measurements  from  the 
previously  laid  out  coordinate  axes  (see  Test  Facilities).  The  stand 
height  was  measured  from  the  floor  level.  The  overhead  lights  were  then 
dimmed  to  start  the  test  run. 

The  observer  sat  in  the  pilot's  seat  with  his  forehead  positioned 
against  the  brow  pad.  He  searched  for  light  images  on  the  canopy  surfaces. 
The  observer  located  three  types  of  light  Images.  The  first  was  a  direct 
view  of  the  light  source  through  the  canopy.  The  second  type  was  a  primary 
order  reflection  in  which  the  light  image  was  reflected  once  from  the 
surface  before  reaching  the  observer's  eye.  The  third  type  was  higher 
order  reflections  in  which  the  light  image  was  reflected  more  than  once 
before  reaching  the  observer's  eye.  The  higher  order  reflections  were 
distinguished  from  the  primary  orders  by  their  movement  when  other  canopy 
windows  were  disturbed. 

The  type  of  data  collected  at  the  end  of  each  test  run  was  determined 
by  the  particular  test  phase.  For  the  first  phase,  the  observer  focused 
the  telephotometer  and  centered  it  on  the  light  image  while  the  assistant 
recorded  the  luminance  value.  This  was  repeated  for  each  of  the  located 
light  images  in  turn. 

For  the  second  phase,  the  observer  used  a  metal  pointer  to  pinpoint 
the  location  of  the  light  image  on  the  canopy  surface.  The  assistant 
marked  the  position  with  a  small  tab  of  adhesive  tape.  He  then  measured 
and  recorded  the  distances  between  the  image  point  and  a  convenient  canopy 
frame  vertex.  The  distances  were  measured  along  the  coordinate  axes  of  the 
mockup.  This  was  repeated  for  each  of  the  located  images  in  turn. 
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RESULTS 


The  reduction  of  the  data  for  the  luminance  and  coordinates  of  the 
test  light  images,  computation  of  a  matched  sample  by  computer  program 
using  the  mathematical  model  for  light  reflections,  and  the  statistical 
analyses  of  the  matched  sample  of  measured  and  computed  values  are 
discussed  below. 

1.  Data  Reduction 

a.  The  transmittance  values  of  the  light  source  were  calculated 
for  the  image  points  of  the  first  test  phase  data.  The  transmittance 
values  equal  the  luminance  readings  of  the  light  images  divided  by  the 
luminance  reading  of  the  light  source. 

b.  The  three-dimensional  coordinates  of  the  imag  surface 
points  in  the  mockup  coordinate  system  were  calculated  for  the  second 
phase  test  data.  The  coordinates  of  the  image  points  equal  the  sum  of  the 
coordinates  of  the  appropriate  vertex  and  the  corresponding  displacement 
distances  which  were  measured  for  the  image  points.  The  coordinates  of  the 
frame  vertices  were  measured  separately  using  a  surveyor's  level  (see 
Apparatus). 


2.  Computation  of  a  Matching  Sample 

a.  A  computer  program  (see  Appendix  C)  was  written  to  compute 
the  transmittance  values  and  coordinates  of  the  images  on  the  canopy 
surface  of  a  light  source.  The  program  was  developed  using  ray  tracing 
concepts  (see  Appendix  B)  and  is  applicable  to  direct  view  and  primary 
order  reflections. 

b.  The  input  data  to  the  program  are  (1)  the  canopy  frame  data, 
(2)  the  location  of  the  mockup  on  the  bay-area  floor,  and  (3)  the  location 
of  the  test  lights  on  the  floor.  The  coordinates  of  the  vertices  of  the 
frame  and  canopy  surfaces  are  program  input  data.  The  vertice  coordinates 
were  measured  after  the  test  using  a  surveyor's  level  (see  Apparatus).  The 
coordinates  of  the  origin  of  the  mockup' s  coordinate  system,  measured  in 
the  rectangular  coordinate  system  on  the  bay-area  floor  (see  Methods),  are 
program  input  data.  The  orientation  of  the  mockup' s  coordinate  system,  as 
measured  from  the  floor  system,  is  also  input  to  the  program.  Finally,  the 
floor  system  coordinates  of  the  light  source  positions  in  the  test  are 
read  by  the  program. 

c.  The  program  computes  the  transmittance  value  and  the 
coordinates  of  the  direct  view  and  reflected  image  points  on  the  canopy 
surfaces  for  each  light  source  position.  The  computations  are  done  for 
both  eyes  of  the  observer,  one  to  the  right  and  the  other  to  the  left  of 
the  pilot's  nominal  eye  position.  The  image  point  coordinates  are  computed 
in  the  coordinate  system  of  the  mockup.  The  printout  of  computed  image 
point  transmittance  and  coordinates  for  each  of  the  light  source  test 
positions  is  a  matching  sample  to  the  values  measured  in  the  mockup  test 
(see  Methods). 
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d.  The  measured  and  computed  values  for  the  image-point 
transmittance  and  coordinates  are  listed  in  Table  1A  for  the  light  source 
test  positions  (see  Appendix  A).  The  table  shows  that  the  coordinates 
computed  for  the  right  eye  vision  are  slightly  closer  in  value  to  the 
measured  values.  This  is  in  agreement  with  the  fact  that  the  observer  is 
right  eye  dominant.  The  right  eye  set  of  data  was  chosen  as  the  matching 
sample. 

e.  The  table  lists  two  computed  image  points  that  were  not 
observed  in  the  test  and  two  observed  values  that  were  not  computed.  A 
physical  check  of  the  mockup  and  the  light  source  locations  for  these 
unmatched  light  images  showed  that  they  were  caused  by  slight 
discrepancies  between  the  mockup  and  the  frame  input  data.  These 
discrepancies  are  listed  in  the  next  paragraph.  The  unmatched  images  were 
dropped  from  the  samples  leaving  two  balanced  and  matching  samples. 

f.  The  unmatched  images  dropped  from  the  samples  and  the 
reasons  for  their  occurrances  are  as  follows: 

(1)  A  computed  image  point  that  was  in  actuality  blocked 
from  the  observer's  view  by  the  brow  pad  which  was  not  included  in  the 
frame  input  data. 

(2)  A  computed  image  point  that  was  blocked  from  view  by  two 
overlapping  canopy  frame  edges  that  were  represented  as  slightly 
separated. 


(3)  An  observed  image  point  which  occurred  at  the  edge  of 
the  canopy  surface  and  which  showed  a  portion  of  the  1-1/2-inch  diameter 
test  bulb,  but  which  was  not  computed  for  the  point  light  source  assumed 
in  the  model. 


(A)  An  observed  image  point  which  did  not  occur  in  the 
computations  because  the  mockup  nose  was  assumed  to  be  2  inches  longer 
than  it  actually  was,  and  the  light  source  was  placed  near  the  nose. 

3.  Statistical  Analyses 

The  difference  score  for  the  transmittance  and  the  three 
coordinates  (x,y,z)  were  treated  in  a  students'  t-test  for  matched  samples 
(Welkowitz  [6],  p.142-143).  The  difference  scores  were  calculated  by 
subtracting  the  measured  values  from  the  computed  values  for  the  matching 
image  points  of  a  given  light  source  position  (see  Appendix  A).  The  main 
difference  scores  (D),  the  standard  deviation  ( Sj-j)  of  the  difference 
scores  and  the  corresponding  t-values  (tp)  are  listed  in  Table  1.  The 
t-values  were  calculated  by  dividing  the  mean  difference  scores  by  the 
standard  error  of  the  mean.  The  null  hypothesis  is  that  the  mean  of  the 
difference  scores  is  zero.  The  degrees  of  freedom  (df)  for  the  number  of 
sample  pairs  are  listed  along  with  the  corresponding  critical  t-values 
(tc)  at  the  .05  level  of  significance  for  a  two-tailed  test  (6,  Table  c, 
p.  257).  The  dimension  of  the  transmission  value  is  in  fractional  parts 
and  that  of  the  coordinate  value  is  in  Inches. 
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TABLE  1 


Student's  t-Test  for  Matched  Samples 


Variable 

D 

SD 

lD 

df 

transmittance 

.0112 

.0471 

.8257 

11 

2.228 

x-coordinate 

-.6378 

.9769 

-2.9197 

19 

2.093 

y-coordinate 

-.6940 

1.4163 

-2.1914 

19 

2.093 

z-coordlnate 

1.6190 

1.1491 

6.3009 

19 

2.093 

*.05  level  of 

significance 

for  a  two- 

tailed  test. 

A  comparison  of  the  t-values  listed  in  Table  1  for  the  scores,  with 
the  critical  t-values,  shows  that  there  is  not  sufficient  reason  to  reject 
the  hypothesis  that  the  measured  and  computed  transmittance  values  are  the 
same.  However,  the  hypothesis  that  the  measured  and  computed  coordinate 
values  are  not  different  must  be  rejected  at  least  at  the  .05  level  of 
significance  (see  Discussion). 

The  matched  data  for  the  transmittance  and  the  three  coordinates  were 
treated  in  linear  regression  analyses  (6,  p.  152-170).  Table  2  lists  the 
Pearson's  correlation  coefficients  (r)  for  the  matched  data,  the  degrees 
of  freedom  (df)  given  by  df-N-2  where  N  is  the  number  of  matched  data,  and 
the  corresponding  critical  r-value  (rc)  at  the  .01  level  of  significance 
for  a  two  tailed  test  (6,  Table  D,  p.  258).  The  table  also  lists  the  slope 
(b)  and  intercept  (a)  of  the  linear  regression  lines  and  the  corresponding 
standard  errors  of  estimate  (se).  ine  coefficients  of  the  linear 
regression  line  were  computed  using  the  method  of  least  squares. 

The  correlation  coefficients  listed  in  Table  2  are  nearly  equal  to 
unity.  A  comparison  of  the  coefficients  with  the  critical  values  shows 
that  the  coefficients  are  statistically  significant  from  zero  at  the  .01 
level.  The  matched  data  are  linearly  related  in  the  statistical  sense. 


DISCUSSION 

It  was  shown  in  the  results  that  the  mean  difference  scores  listed  in 
Table  1  for  the  x-,  y-,  and  z-coordinates  are  statistically  significant 
from  zero.  The  mean  difference  scores  are  equal  to  the  differences  between 
the  sample  means  for  the  measured  and  computed  coordinates.  The  mean 
differences  may  be  the  errors  for  centering  and  aligning  the  coordinate 


TABLE  2 


Linear  Regression  Analysis  for  the  Match  Data 


Variable 

r 

df 

rc* 

b 

a 

8e 

t  ransmittance 

.9923 

10 

.708 

.9476 

.009 

.0444 

x-coordinate 

.9985 

18 

.561 

.9145 

.3561 

.9228 

y-coordinate 

.9973 

18 

.561 

.9987 

.8385 

1.380 

z-coordinate 

.9897 

18 

.561 

1.0981 

-10.655 

1.149 

*.01  level  of  significance  for  a  two-tailed  test. 


system  of  the  test-area  floor  with  that  of  the  mockup.  The  centering  was 
done  by  plumb  bob  and  the  alignment  by  tape  measure  and  line  of  sight.  The 
height  of  the  mockup' s  top  canopy  surface  above  the  test  area  floor  was 
measured  by  tape  sometime  after  the  test.  The  Intervening  period  was  used 
by  other  personnel  to  measure  the  coordinates  of  the  vertices  of  the 
mockup  frame.  The  mockup  was  Inadvertently  braced  during  the  operation 
thereby  changing  slightly  Its  height  above  the  floor  level. 

The  90  percent  confidence  limits  for  the  slopes  (b)  and  the  intercepts 
(a)  of  the  linear  regression  lines  are  listed  In  Table  3  for  the  test 
variables.  The  table  lists  also  the  standard  estimates  of  error  (se)  for 
the  measured  variables  and  the  corresponding  degrees  of  freedom  (df-N-2). 
It  Is  assumed  that  the  measured  values  for  a  given  computed  value  form  a 
normal  distribution  (1,  p.  197-198).  The  confidence  limits  for  the  slope 
and  Intercept  are  computed  using  the  standard  error  of  the  mean  and  the 
critical  t-value  (.05  level  significance)  for  the  degrees  of  freedom  (1, 
Table  A5 ,  p.  464).  The  confidence  bounds  for  the  standard  errors  of 
estimate  are  computed  using  the  critical  values  (  .05  and  .95  level  of 
significance)  of  the  ratio  of  Chl-squared  to  degrees  of  freedom  (1,  Table 
A6b,  p.  466). 

The  slope  of  the  linear  regression  line  should  be  equal  to  unity  (b-1) 
and  the  intercept  equal  to  zero  (a-0)  when  the  matched  samples  are  equal 
within  the  limits  of  the  measurement  error.  Furthermore,  the  standard 
error  of  measurement  should  be  equal  to  the  measurement  error  for  the 
equipment  and  procedure.  These  values  are  contained  within  the  confidence 
limits  for  the  transmittance  as  shown  In  Table  3. 
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TABLE  1 


Student's  t-Test  for  Matched  Samples 


Variable 

D 

SD 

CD 

df 

t  * 
Lc 

transmittance 

.0112 

.0471 

.8257 

11 

2.228 

x-coordinate 

-.6378 

.9769 

-2.9197 

19 

2.093 

y-coordinate 

-.6940 

1.4163 

-2.1914 

19 

2.093 

z-coordinate 

1.6190 

1.1491 

6.3009 

19 

2.093 

*.05  level  of 

significance 

for  a  two- 

tailed  test 

• 

A  comparison  of  the  t-values  listed  In  Table  1  for  the  scores,  with 
the  critical  t-values,  shows  that  there  is  not  sufficient  reason  to  reject 
the  hypothesis  that  the  measured  and  computed  transmittance  values  are  the 
same.  However,  the  hypothesis  that  the  measured  and  computed  coordinate 
values  are  not  different  must  be  rejected  at  least  at  the  .05  level  of 
significance  (see  Discussion). 

The  matched  data  for  the  transmittance  and  the  three  coordinates  were 
treated  in  linear  regression  analyses  (6,  p.  152-170).  Table  2  lists  the 
Pearson's  correlation  coefficients  (r)  for  the  matched  data,  the  degrees 
of  freedom  (df)  given  by  df=N-2  where  N  is  the  number  of  matched  data,  and 
the  corresponding  critical  r-value  (rc)  at  the  .01  level  of  significance 
for  a  two  tailed  test  (6,  Table  D,  p.  258).  The  table  also  lists  the  slope 
(b)  and  intercept  (a)  of  the  linear  regression  lines  and  the  corresponding 
standard  errors  of  estimate  (se).  The  coefficients  of  the  linear 
regression  line  were  computed  using  the  method  of  least  squares. 

The  correlation  coefficients  listed  in  Table  2  are  nearly  equal  to 
unity.  A  comparison  of  the  coefficients  with  the  critical  values  shows 
that  the  coefficients  are  statistically  significant  from  zero  at  the  .01 
level.  The  matched  data  are  linearly  related  in  the  statistical  sense. 


DISCUSSION 

It  was  shown  in  the  results  that  the  mean  difference  scores  listed  in 
Table  1  for  the  x-,  y-,  and  z-coordinates  are  statistically  significant 
from  zero.  The  mean  difference  scores  are  equal  to  the  differences  between 
the  sample  means  for  the  measured  and  computed  coordinates.  The  mean 
differences  may  be  the  errors  for  centering  and  aligning  the  coordinate 


TABLE  2 


Linear  Regression  Analysis  for  the  Match  Data 


Variable 

r 

df 

rc* 

b 

a 

8e 

t  ransmittance 

.9923 

10 

.708 

.9476 

.009 

.0444 

x-coordlnate 

.9985 

18 

.561 

.9145 

.3561 

.9228 

y-coordinate 

.9973 

18 

.561 

.9987 

.8385 

1.380 

z-coordlnate 

.9897 

18 

.561 

1.0981 

-10.655 

1.149 

*.01  level  of  significance  for  a  two-tailed  test. 


system  of  the  test-area  floor  with  that  of  the  mockup.  The  centering  was 
done  by  plumb  bob  and  the  alignment  by  tape  measure  and  line  of  sight.  The 
height  of  the  mockup* s  top  canopy  surface  above  the  test  area  floor  was 
measured  by  tape  sometime  after  the  test.  The  intervening  period  was  used 
by  other  personnel  to  measure  the  coordinates  of  the  vertices  of  the 
mockup  frame.  The  mockup  was  inadvertently  braced  during  the  operation 
thereby  changing  slightly  its  height  above  the  floor  level. 

The  90  percent  confidence  limits  for  the  slopes  (b)  and  the  intercepts 
(a)  of  the  linear  regression  lines  are  listed  In  Table  3  for  the  test 
variables.  The  table  lists  also  the  standard  estimates  of  error  (se)  for 
the  measured  variables  and  the  corresponding  degrees  of  freedom  (df-N-2). 
It  Is  assumed  that  the  measured  values  for  a  given  computed  value  form  a 
normal  distribution  (1,  p.  197-198).  The  confidence  limits  for  the  slope 
and  Intercept  are  computed  using  the  standard  error  of  the  mean  and  the 
critical  t-value  (.05  level  significance)  for  the  degrees  of  freedom  (1, 
Table  A5 ,  p.  464).  The  confidence  bounds  for  the  standard  errors  of 
estimate  are  computed  using  the  critical  values  (  .05  and  .95  level  of 
significance)  of  the  ratio  of  Chl-squared  to  degrees  of  freedom  (1,  Table 
A6b,  p.  466). 

The  slope  of  the  linear  regression  line  should  be  equal  to  unity  (b"l) 
and  the  intercept  equal  to  zero  (a«0)  when  the  matched  samples  are  equal 
within  the  limits  of  the  measurement  error.  Furthermore,  the  standard 
error  of  measurement  should  be  equal  to  the  measurement  error  for  the 
equipment  and  procedure.  These  values  are  contained  within  the  confidence 
limits  for  the  transmittance  as  shown  in  Table  3. 


TABLE  3 


Confidence  Limits  (90%)  for  the  Slope  (b)  and  the  Intercept  (a) 
of  the  Linear  Regression  Line,  and  the  Standard  Errors  of 

Estimate  (sc) 


Slope  (b) 

Intercept  (a) 

Error 

(se) 

Variable 

df 

max 

min 

max 

min 

max 

min 

transmittance 

10 

1.0119 

.8834 

.0322 

-.0143 

.1252 

.0328 

x-coordinate 

18 

.9403 

.8887 

.7139 

-.0017 

1.8567 

.7695 

y-coordinate 

18 

1.0272 

.9702 

1.3736 

.3034 

2.7767 

1.091 

z-coordinate 

18 

1.1603 

1.0359 

-10.2098 

-11.1008 

2.3119 

.9084 

The  confidence  limits  in  Table  3  for  the  coordinates  do  not  in  all 
cases  contain  the  expected  values  noted  above  for  the  regression  line 
slopes  and  intercepts.  The  confidence  limits  for  the  slopes  of  the  x  and  z 
coordinates  are  close  to,  but  do  not  contain  the  expected  value  (b=l). 
Similarly,  the  limits  for  the  y  coordinate's  intercept  are  close  to  the 
expected  value  (a-0)  but  do  not  contain  it.  This  discrepancy  is  highly 
noticeable  In  the  case  of  the  z  coordinate's  intercept  which  is  on  the 
order  of  -10  Inches.  However,  the  values  of  the  z  coordinate  range  from  80 
to  100  inches.  The  line  intercept  must  correct  for  an  additional  increment 
in  value  due  to  the  slope  (or  the  order  of  0.1x90*9  inches),  as  well  as 
the  1.62  inches  mean  difference  between  the  computed  and  measured  samples. 
A  slight  tilt  of  the  mockup  on  the  bay  area  floor,  not  accounted  for  in 
the  alignment  measurements,  is  a  reasonable  explanation  for  the  slope 
values  of  the  x  and  z  coordinate  samples.  The  observer  would  dismount  from 
the  mockup  during  the  second  test  phase  to  assist  in  coordinate 
measurements.  The  resulting  weight  shift  in  the  mockup  between  the  times 
of  observing  and  measurement  could  be  the  cause  of  a  tilt  in  the  mockup 
f  rame . 

Similarly,  the  lower  limits  on  the  standard  errors  of  estimate  are 
larger  (on  the  order  of  1  inch)  than  would  be  expected  for  the  equipment 
and  procedure  used  to  measure  the  positions  of  the  light  image  (see 
Apparatus).  However,  an  in-depth  error  analysis  would  include  the  effects 
of  errors  in  locating  the  light  source  on  the  bay  area  floor  and  in 
pinpointing  the  position  of  the  image  of  the  l?5-inch  diameter  light  bulb. 
These  additional  factors  may  explain  the  otherwise  larger-than-expected 
standard  errors  of  estimate. 
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RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


The  following  Is  recommended: 

1.  The  model  should  be  tested  using  cylindrically  shaped  panels  as 
well  as  the  planar  panels  used  in  this  test.  The  use  of  the  surveyor’s 
level  to  measure  coordinate  values  should  be  extended  to  measuring  the 
positions  of  the  light  source  and  the  surface  image  points  as  well  as  the 
mockup  frame  vertices. 

2.  The  locations  on  the  canopy  surface  of  the  primary  reflections 
of  the  Instrument  panel  lights  can  be  computed  using  the  computer  model. 
This  area  of  application  should  be  Investigated  further. 


CONCLUSION 

The  US  Army  Human  Engineering  Laboratory  (USAHEL)  has  experimentally 
verified  a  computer  model  for  the  Internal  light  reflections  on  the 
transparent  surfaces  of  helicopter  canopies.  The  model  was  verified  using 
a  mockup  of  the  Model  209  AH-1S  Helicopter  with  the  flat  plate  canopy 
design.  The  transmittance  values  and  the  coordinates  of  the  light  images 
on  the  canopy  surfaces  were  measured  for  a  wide  range  of  light  source 
positions.  A  matched  sample  was  then  computed  for  the  source  positions 
using  the  computer  model. 

The  matched  values  of  the  measured  and  computed  samples  show  a  high 
correlation.  The  Pearson's  correlation  coefficients  for  a  linear 
regression  analyses  are  greater  than  0.98  and  the  coefficients  are 
significantly  different  from  zero  at  the  .01  level.  However,  the  mean 
differences  between  the  measured  and  computed  coordinate  samples  are  on 
the  order  of  0.6  Inches  to  1.62  Inches.  Furthermore,  the  slopes  and 
intercepts  of  the  linear  regression  lines  are  not  the  expected  values  for 
all  of  the  coordinate  samples.  These  differences  are  presumably  due  to  the 
Inaccuracies  In  the  alignment  of  the  mockup  with  the  coordinate  system  on 
the  test-area  floor. 
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TEST  DATA  AND  STATISTICAL  ANAYSES 


The  measured  and  computed  values  of  the  light  image  transmittance  and 
coordinates  are  listed  in  Table  1A  for  the  test  positions  of  the  light 
source.  The  table  lists  the  transmittance  and  coordinates  of  the  matched 
direct  view  and  primary  reflection  data  for  the  first  phase.  The  second 
phase  data  Includes  the  coordinates  of  the  matched  direct  view  and  primary 
reflection  data,  as  well  as  the  unmatched  reflection  data. 

The  light  image  coordinates  are  in  the  coordinate  system  (xm,  ym,  zm) 
of  the  mockup,  which  was  aligned  as  closely  as  possible  to  the  design 
specifications  for  the  actual  aircraft.  The  coordinates  of  the  light 
source  positions  are  in  the  coordinate  system  (xf,  yf  ,  z^ )  laid  out  on  the 
test  area  floor.  The  origin  of  the  mockup' s  system  was  located  in  the 
floor  system  as  follows:  Xf  ■  o,  y^  -  24.75,  and  Zf  •  -30.586  inches.  The 
z  axis  of  the  mockup  was  directed  parallel  to  that  of  the  floor  system, 
while  the  y  axis  was  directed  in  the  parallel  but  opposite  direction  of 
that  of  the  floor.  A  light  source  located  at  the  coordinates  Xf,  yf ,  zf  in 
the  floor  system  was  located  at  the  mockup  coordinates:  xm  •  -X{  ,  ym  ■  -y{ 
+  24.75  and  zm «  zf  +  30.586  inches. 

The  attached  computer  program  MSTT  was  used  to  compare  the  matched 
data  in  a  student's  t-test.  The  data  were  treated  in  a  linear  regression 
analysis  using  the  program  LSA  attached  below. 
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TABLE  1A 


Transaictance  (T)  and  Position  Coordinates  (x,y,z)  of  the  Direct  View  (D) 
and  Primary  (P)  Reflected  Images  for  the  Light  Source  Test  Positions 


Light  Source  _ Computed _  _ Measured 


X 

I. 

y  z  Type  xc 

First  Test  Phase  Data 

z 

c 

T_c 

z_m 

T 

m 

60. 

-52. 

12. 

P  -7.65 

130.14 

105.97 

.0734 

.0901 

36. 

-40.4 

12. 

P  -5.96 

128.49 

105.97 

.0543 

.0502 

36. 

7.56 

36. 

D  -14.77 

87.06 

83.12 

.7235 

.7907 

60. 

43.57 

24. 

D  -15.29 

96.83 

84.5 

.7672 

.7326 

P  15.34 

111.74 

91.86 

.1019 

.0858 

48. 

-88.4 

48. 

D  -16.75 

129.62 

89.31 

.8693 

.9002 

P  15.14 

133.53 

95.36 

.0658 

.1089 

72. 

-64.4 

48. 

D  -16.38 

128.17 

91.76 

.8672 

.7698 

P  15.01 

131.83 

96.2 

.0663 

.0494 

72. 

-4.43 

48. 

D  -16.45 

115.07 

91.74 

.8239 

.7531 

P  14.88 

121.62 

96.28 

.076 

.0673 

12. 

31.56 

24. 

P  -2.84 

114.26 

105.98 

.1347 

.0908 

II. 

Second  Test 

Phase  Data 

• 

00 

-100.44 

36. 

D  -17.45 

133.78 

84.5 

- 

16. 

133.69 

83.17 

P  -9.31 

136.09 

105.96 

Note  (1) 

60. 

-88.44 

36. 

D  -17.11 

131.31 

86.85 

- 

15.75 

131.31 

85.16 

P 

13.78 

136.91 

91.81 

Note  (2) 

36. 

-76.44  12. 

P  -4.24 

133.61 

105.96 

-2.72 

133.22 

105.97 
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TABLE  1A  (Continued) 


Transmittance  (T)  and  Position  Coordinates  (x,y,*)  of  the  Direct  View  (D) 
and  Primary  (P)  Reflected  Images  for  the  Light  Source  Test  Positions 


Light  Source 


Computed 


Measured 


X 

y 

z 

Type  xc 

L£ 

!c 

T_c 

>m 

z_ 

m 

72. 

-40.44 

48. 

D 

-16.48 

121.30 

91.38 

-15.25 

121.94 

88.91 

P 

14.97 

128.50 

96.20 

14.28 

130.66 

93.31 

84. 

-40.44 

24 

D 

-17.18 

123.19 

86.67 

-15.75 

123.69 

84.66 

P 

15.31 

129.80 

93.67 

14.66 

132.16 

91 . 31 

48. 

-28.44 

24. 

D 

-18.40 

104.99 

79.24 

-16.12 

107.94 

78.16 

P 

-7.29 

124.43 

105.97 

-5.72 

123.97 

105.97 

36. 

-4.43 

36. 

D 

-15.11 

91.15 

82.85 

-15.5 

90.69 

80.41 

P 

-6.75 

115.9 

105.98 

-5.22 

115.47 

105.97 

P 

15.34 

115.75 

92.24 

14.782 117.16 

89.44 

72. 

7.57 

48. 

D 

-16.54 

109.71 

91.36 

-15.5 

111.94 

88.91 

P 

14.89 

121.78 

96.18 

14.532  123.16 

93.31 

12. 

19.57 

12. 

P 

-5.34 

73.22 

92.81 

36. 

31.57 

48. 

D 

-13.53 

81.72 

88.8 

-12.78 

81.97 

87.69 

P 

14.01 

95.89 

94.2 

14.782 

99.66 

92.31 

60. 

31.57 

24. 

P 

15.35 

116.45 

92.24 

14.782 

117.16 

90.31 

12. 

43.57 

48. 

D 

-4.8b 

67.87 

87.83 

-3.375 

68.94 

88.75 

24. 

43.57 

36. 

P 

13.61 

82.38 

88.62 

13.78 

83.41 

87.81 

60. 

43.57 

48. 

D 

-14.46 

98.97 

91.14 

-15. 

96.44 

88.41 

P 

14.86 

113.75 

95.79  Note(4) 

NOTES  -  An  examination  of  the  mockup  showed  the  following  reasons  for  the 
occurrence  of  unmatched  data: 

(1)  The  computed  reflected  image  was  actually  blocked  from  the  observer's 
view  by  the  brow  pad  which  was  not  included  in  the  mockup  frame  vertice  data. 

(2)  The  observed  image  of  the  IV inch  diameter  bulb  appeared  at  the 
window  edge  and  was  not  conputed  for  the  point  source  used  in  the  model . 

(3)  The  mockup  nose  was  assumed  to  be  2  inches  longer  in  the  model  than 
it  actually  was,  and  the  light  source  was  placed  just  under  the  nose. 

(4)  The  computed  image  was  actually  blocked  from  view  by  two  overlapping 
frame  edges,  but  they  were  represented  as  slightly  separated  in  the  model. 


(Concluded) 


SMYTH*  STMFZ . 

ACCOUNT*  HE*** • 

FTN  . 
l  GO. 

PROGRAM  MSTT( INPUT .OUT PUT*  TAPE 2 -INPUT* TAPE  1 -OUT PUT! 

C  T-TE ST  FOR  MATCHED  SAMPLES 
C  XM.  RANDOM  VARIA8LE  MEASUREO  BY  TESTING 
C  XC»  COMPUTED  VALUES  FUR  TEST  MEASUREMENT  CONDITIONS 
C  N*  NUMBER  OF  PAIRED  OBSERVATIONS 

DIMENSION  XMd00)*XC<100)*FLd00)*D(l00) 

RE  AC ( 2  *  1000 )  NR 
DO  100  IR ■ 1  *  NR 

C  RE AO  IN  TEST  AND  COMPUTEO  VAR  IB L  c  VALUES 
RE AD (2, 999) (Fit  I ) » I -1* 40 ) 

999  FORMAT ( AO A2 ) 

WRITE ( 1*999) <FL< I)» 1-1,40) 

READ(2.1000)N 

1000  F0RMAT(2X,I4) 

XN-N 

READ(2»1001)(XM(I)*I«1»N) 

1001  F0RMAT(2X,6(F10.4*2X) ) 

RE AO (2*1001 ) t  XC ( I ) *  1-1*  N) 

WRITE(  1  *  101 1 ) N 

1011  F0RMAT(2X»'TEST  DATA1 /2X» 'NO.  MEASUREMENTS- •» 14 /2X» • TEST  MEASUREME 
QNT  S  •  ) 

WRITE (1*1012) l XM1I) *I«1,N) 

1012  FORMAT <6(2X»F10.4) I 
WRITE(1» 1013) (XC(l). 1-1. N) 

1013  F0RMAK2X,  'COMPUTED  V  ALUE  S  •  /  6(  2X,  F  10 . 4  )  ) 

DO  5  I-l.N 

0 ( I ) -XC ( I  )“XM ( I ) 

5  CONTINUE 

WRITE11.1015) (0(1)* I-l.N) 

1015  FORMAT (2X, "DIFFERENCE S"/6(2X,F 10.4) ) 

SD-O. 

SD2-C. 

DO  10  I - 1 »  N 
SD-SD+D ( I ) 

S02«SD2+D(I)*D(I) 

10  CONTINUE 
SMD-SO/XN 

ST D- SORT ( CS02-SD*SD/XN ) / ( XN-1. ) ) 

TV-SD/S0RTC(XN*SD2-SD*S0)/(XN~1. ) ) 

WRITE (1*1014) SMD.STD* TV 

1014  FORMAT ( 2X ,*Mt AN-DIF FE R ENC E", 2 X , F 10. 4, 2X* "STANDARD  DEVIATION", 
Q2X*F10.4/2X,"T-TEST  V ALUE ", 2X, F 1 0 . 4 ) 

100  CONTINUE 
STOP 
ENO 
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SHY  T  H»  ST  Hf  Z  . 

ACC  OUN T#  HE  *  **  . 

FTN. 

160. 

PROGRAM  LSA( INPUT  *  OUT  PUT  » T A P t  2- INPUT* T AP c 1» OU'TP UT I 
C  LINEAR  REGRESSION  ANALYSIS  BY  LEAST  SQUARES 
C  XM,  KANCCM  VARIABLE  ME  ASURtD  BY  TESTING 
C  XCt  CCMPUT  £  0  VALUES  FUR  TEST  MEASUREMENT  CONDITIONS 
C  N,  NUMBE*  OF  PAIRED  OBSERVATIONS 

DIMENSION  XMI 100 >/ XC t 1 00 >,FL 11001 
PE  AO ( 2.  1 000 ) NR 
00  100  I P  *  1 »  N R 

C  READ  IN  TtST  ANO  COMPUTED  VARI3LE  VALUES 
P E A D ( 2 , 9 9 9 ) I F L (II . i - 1 . 4 0  ) 

999  f  OF1*/ T  <  4GA2  ) 

WRITE  (l»999)(Fl<I  1.  1*1.  4J) 

KEAC(  2  ,  10001N 

1000  FORMAT ( ?X . I  4 ) 

X  N  *  N 

RE  AC (2. LG01  I  ( XM  (  I  ),  1*1, N) 

1001  FORMAT (2X»6(F10,4»2X)  ) 

RE  ADI ?. 1001 ) (XC (I), I-l.N) 

WRITE(1,1011)N 

1011  EuRMAT<?X,'TeST  DATA' /2X, 'NO.  ME  A  SURE  ME  NT S * » » I  4 / 2 X. » TE S T  MEASUREME 
QMS'  ) 

WR  ITE  (1 , 1C12MXMI  I »  ,  I-l.N) 

1012  FORMAT <6<2X, FID. 4>) 

WR  ITE  ( l»  1013HXC  <  I  > .  I-l.N) 

1013  F0RMAK2X, 'COMPUTED  V ALUE S ' / 6 < 2X> F 10 . 4 ) > 

C  READ  fuR  C0NHLEN1  INTERVAL  CALCULATIONS  THE  TYPE  1  ERROR  VALUE  AND 
C  CORRESPONDING  TWO-SIDED  STUDENT'S  T-VALUE  AT  N-2  DEGREES  CE  FREEDOM") 
RE  AD ( 2. 1010) AT, TV 
1010  FORMAT (2X,3(F10.4,2X)  ) 

C  READ  FOR  TEST  VARIANCE  CALCULATIONS  TYPE  1  ERROR  LEVEL  ANO  CCRkESPDNDI 
C  RATIO  OF  CHI-SQUAREO  TO  DtCREcS  OF  FRELDOM 
P  E  A  0 ( 2  » 10 10  I  AC. CPU. CF  l 

C  COMPUTE  LINEAR  CORRELATION  COEFFICIENT  FOR  PAW  MEASURED  AND  COMPUTED  V 
SMC-O. 

SM-C. 

$C*C. 

SC  2*C . 

SM2-C. 

00  1C  I-l.N 
SM-SM+XMl  1 ) 

SC- SC RXC (  I ) 

SMC«SMCmXM(  I  )*xc (  I  ) 

SM2»SM2*XM( 1  )  * X M (  I  ) 

SC2-SC2«-XCm*XC<  l ) 

1C  CONTINUE 

SSC-SQRT ( (SC2-SC*SC  A  XN ) / ( XN-1,  )  ) 

RXY2-I XN*SMC-SC*SM) **2 

RXY2-RXY2/(<  <CN*SC  2-SC*SC  >  *1  XN*  SM2-SH*SM  >  ) 

RXY*SQRT( RXY2  ) 

WRITE (1.1C02)RXY»RXY2 

1002  F0RMATI2X, 'CORRELATION  C  QE  F  F  1C  I  ENT  »  '  » E  7 . 4 , 2  X ,  '  S  QU  AP  fc  C  COEFFICIENT- 
Q*»F 10.41 

C  COMPUTE  LINEAR  REGRESSION  EQUATION  COEFFICIENTS  FOR  ESTIMATING  MEASURE 
C  XH1ESTIMATE0) »  8*XC  +  A 
AVM-SM/XN 
AVC-SCXXN 
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B»RXY*<XNASM2-SM*SM)/<XN*SC2-SC*SCI 

A«AVM-8*A VC 

WRITE<1,1003>R#A 

1003  FGRMATI2X, 'REGRESSION  EOUATION  COE  F  E  IC  I  E  N 1  S  *  B «  •  » F 10  .A , 2X»  •  A*  •  , F  10 
O.A  ) 

C  COMPUTE  STANDARD  tRRDR  OF  tSTIMATE 
SE-O. 

00  20  I  •  1  *  N 

SE«SEF(Xhm-BPXC  ( I)-A1**2 
20  CONTINUE 

SE-SORT(SE/XNI 
WR ITE  < 1, 1C0A  )  SE 

100 A  FORMAT { 2  X  , • STANDARD  IRRUR  OF  E S T 1 M A T E ■ • *  El 0  .  A  ) 

C  COMPOTE  CONFIDENT  1NTERVAI  FOR  REGRESSION  EOUATION  COEFFICIENTS 
RGA.TV*SE /SORT ( XN) 

RGB«TVASF /(SORT ( XN-1. I*SSC) 

BL • [-RGB 
BU-o+ROB 
AL ■ A-RGA 
AU»  A  +  KG A 
wRITr < 1, 1005IA1 

1005  FOR  MAT l ?X * '  CONF ID ENT  INTERVAL  TYPE  1  ERROR  LEVEL*  A«»*F10.A) 

WRITE! 1* 1006) BL, BU 

ICOt  FORMAT! 2X, 'REGRESSION  EOUATION  CUE F F  IC 1 E NT • / 2 X,  •  L Uwt R  LIMIT,  BL  »  ' , 
OFIC.A, •  UPPtR  LIMIT,  BU-',F10.A) 

WRI1E(1,1L07)AL,AU 

1007  FORMAT! 2X, ' LOWER  LIMIT,  AL*',FlU.A,'  UPPER  LIMIT,  AU-»,F10.A) 

C  COMPUTE  ACCEPTIBLE  TE S T-VAR I ANC t  DUE  TO  TtST  PROCttOUKE  AND  INSTRUMENT 
STU-SE/SORTICFUI 
S  T  L • SE / SUR  T (CFL  » 

WRITE ( 1,1C09)AC,STU,STL 

1 00 A  F0PMATI2X"TEST  VARIANCE  TYPE  1  EPRLR  LEVEL,  A ■ " , F 1 0  .  A, 2 X / 

Q2X,” TEST-VARlANCf  UPPER  LIMIT,  STU » F 1 0  .  A  /  2 X  , 

0  "L  OWL  R  LIMIT,  S  T  L  - ,  F  1  0  .  A  I 
100  CONTINUE 
ST  CP 
END 
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APPENDIX  B 

PRIMARY  ORDER  REFLECTIONS 


PRIMARY  ORDER  REFLECTIONS 


The  computation  of  the  transmittance  and  position  coordinates  of  the 
direct  view  and  primary  order  reflection  images  are  discussed  in  this 
section.  The  computations  for  a  direct  view  image  are  discussed  first.  A 
primary  order  reflection  occurs  when  the  light  ray  from  the  light  source 
is  reflected  once  before  reading  the  observer's  eye.  We  consider  next  the 
computations  for  the  primary  order  reflections  from  a  planar  surface,  then 
that  from  a  cylindrical  surface. 


DIRECT  VIEW  IMAGES 

The  computational  procedures  for  direct  view  image  points  are  derived 
from  established  techniques.  The  position  coordinates  are  first  computed 
for  the  point  where  the  straight  line  ray  from  the  light  source  to  the 
pilot's  eye  position  intercepts  a  window  surface.  The  directional  cosines 
are  computed  for  the  ray  from  the  end  position  coordinates  and  the  length. 
The  Interception  point  is  computed  from  the  light  source  position,  the 
directional  cosines  of  the  ray  and  the  parameters  of  the  surface  using 
established  techniques  for  a  planar  surface  (Smyth,  June,  1977,  App  A)  and 
a  cylindrical  surface  (Smyth,  July  1977,  App  A).  The  interception  point  is 
tested  against  the  window  edges  to  ensure  enclosure.  The  light 
transmittance  is  computed  next  from  the  angle  of  incidence  between  the  ray 
and  the  surface  normal  (Smyth,  June  1977,  App  A). 


PRIMARY  ORDER  REFLECTIONS  FROM  A  PLANAR  SURFACE 

The  coordinates  of  the  reflection  point  on  the  window  surface  are  a 
function  of  the  light  source  position,  the  pilot's  eye  position,  the 
surface  normal  and  the  position  of  a  window  vertex.  We  consider  the 
geometry  (see  Figure  IB)  formed  by  the  positions  of  the  pilot's  eye  (P0  ■ 
Xo.Yo»*o>i  the  light  source  (P|  -  x(,ya,z()  and  the  reflection  point  (Pr  - 
xr,yT,zr).  The  plane  containing  these  three  points  contains  also  the 
surface  normal  (an,bn,cn)  by  the  laws  of  reflection.  Furthermore,  the 
angle  between  the  ray  P*  Pr  and  the  surface  normal  is  equal  and  opposite  to 
that  between  the  ray  P0Pr  and  the  surface  normal  by  the  same  law. 

The  straight  line,  P,  PQ ,  from  the  light  source  to  the  pilot's  position 
is  specified  by  its  starting  position  (P(  :  x,  , y,  the  directional 

cosines  (Sp,  ,h0l  ,coa )  and  its  length  (Roa).  The  last  two  parameter  sets 
are  easily  computed  from  the  known  end  positions.  We  denote  the  point 
where  the  line  P» P0  Intersects  the  extension  of  the  surface  normal  by  Pj(X|, 
y)(S|).  The  extension  separates  the  triangle  P(  Pr  P0  into  two  smaller 
triangles  with  sides  a  and  b  such  that  a  +  b  ■  R0(.  We  find  it  useful  to 
derive  the  coordinates  of  the  point  (Pj)  as  an  intermediate  step  in  the 
derivation. 


The  coordinates  of  the  point,  Pj ,  are  determined  from  the  ratio  of  the 
length  of  the  line  PgPr  to  that  of  the  line  P,  Pr  ;  i.e.,  Q  ■  RB/Rr.  The 
ratio  is  computed  using: 

1.  the  cosine  of  the  angle  between  the  line  PrPe  and  the  surface 

normal* 

2.  that  of  the  line  PrP8  and  the  surface  normal,  and 

3.  that  of  the  line  PrP0  and  the  normal. 

The  first  term  is  equal  to  zero  while  the  second  and  third  are  equal  to 
each  other.  The  results  are: 

an^xe  “  xs)  +  bn  ^e  “  ys^  +  cn^2e  ~  *s^  (1) 

Q  *  — . — . * -  ■  -  --I 

an(xe  '  xo>  +  bn  <ye  ~  ? o )  +  cn<ze  '  zo> 

The  length  (a)  of  the  line  P8Pe  is  determined  by  the  law  of  sines  for  the 
two  triangles  and  the  expression  Rqs  ■  a  +  b, 

•  “  QRqs  (2) 

T+q" 


The  coordinates  of  the  point,  Pt,  are  given  by: 


xl 

“  X  8 

+ 

a  • 

aos» 

-  y8 

+ 

a  • 

bOS* 

(3) 

*1 

“  *8 

+ 

a  • 

COB* 

The  coordinates  of  the  primary  reflection  point  Pr  are  determined  by 
the  point  Pj ,  the  directional  cosines  of  the  surface  normal,  and  the 
distance,  Rjr,  from  the  point,  Pj,  to  the  surface  in  the  direction  of  the 
surface  normal.  This  distance  may  be  computed  by  established  techniques 
(Smyth,  June  1977,  App  A). 

Rir  •  «n(xi-xe)  +  bn  <yj-ye)  +  cn  (z|-ze>*  (4) 

The  coordinates  of  the  reflection  point  (Pr)  are, 

xr  -  -  *nRir» 

Yr  -  Yi  -  Mir.  (5) 

zr  -  zj  -  cnRir» 

The  point,  Pr,  is  a  reflection  point  for  the  light  source  if  the  value  of 
Q  is  given  by  equation  (1)  is  greater  or  equal  to  zero;  i.e.,  Q>o,  for 
then  the  light  source  and  observer  are  on  the  same  side  of  the  panel 
surface. 


The  reflection  point  (Pr)  is  tested  against  the  window  edges  to  ensure 
enclosure.  The  light  reflectance  is  computed  from  the  angle  of  incidence 
between  the  ray  and  the  surface  normal  (Smyth,  June  1977,  App  A). 


PRIMARY  ORDER  REFLECTIONS  FROM  A  CYLINDRICAL  SURFACE 

The  computation  of  the  primary  reflection  points  for  a  cylindrical 
surface  will  be  derived  first  for  a  simple  case  and  then  extended  to  all 
cases  in  general.  We  consider  a  cylindrical  surface  with  the  following 
parameters  (see  Figure  2B):  (1)  the  origin  for  the  cylindrical  axis  is 
located  at  the  origin  of  the  coordinate  system,  (2)  the  cylindrical  axis 
is  directed  along  the  y  axis  of  the  coordinate  system,  and  (3)  the 
cylindrical  radius  is  equal  to  unity.  The  analysis  is  similar  to  that 
developed  above  for  the  planar  surface. 

Consider  the  geometry  formed  by  the  positions  of  the  light  source,  the 
reflection  point  and  the  pilot's  eye.  The  straight  line  from  the  light 
source  to  the  eye  is  specified  by  its  starting  position,  the  directional 
cosines  (aos,bos,cos)  an<*  t^e  *-ine  leogib  (Ros)*  7116  8urface  normal  (an,bn, 
cn)  is  contained  within  the  plane  defined  by  the  three  points  Ps,Pr 
and  P0 .  We  let  the  point  where  the  extension  of  the  surface  normal 
intersects  the  line  P8P0  be  denoted  by  Pt  and  its  coordinates  by  Xj,  y1 , 
and  zj.  The  line  separates  the  triangle  Psprpo  into  two  smaller  triangles 
and  the  PSP0  into  two  line  segments  of  lengths  a  and  b.  The  length  (a)  of 
the  line  PsPj  is  determined  by  the  law  of  sines  for  the  two  triangles  and 
the  expression  Ros  *  a4b.  The  length  (a)  is  given  by  the  equation  (2), 
where  Q  is  the  ratio  of  the  length  of  the  line  PsPr  to  that  of  the  line  Pr 
PQ;  i.e.,  Q  ■  Rs/R0.  The  coordinates  of  the  point  Pj  are  given  by 
equations  (3). 

The  surface  normal  is  orthogonal  to  the  cylindrical  axis,  and  the 
distance  (R)  along  the  surface  normal  from  the  point,  P|  ,  to  the  point 
(o*Yj»o)  on  the  cylindrical  axis  is  given  by, 

R  -  ((x8  +  Qx0)2  +  (zs  +  Qz0)2)VS  /(I  +  Q).  (6) 

The  directional  cosines  of  the  surface  normal  at  the  reflection  point  may 
now  be  computed  as, 


an  -  (xs  +  Qx0)/R, 

bn  -  0,  (7) 

cn  "  <zs  +  Qzo)/R* 


Finally,  since  the  cylinder  radius  is  unity,  the  coordinates  of  the 
reflection  point  (Pr)  are: 


xr  '  “an» 

yr  -ys+a  •  bo8.  (8) 


The  Q  factor  used  In  equations  (2)  and  (3)  is  determined  by  setting 
the  cosine  of  the  angle  between  the  line  PrP8  and  the  surface  normal, 
equal  to  that  of  the  angle  between  the  line  PrP0  and  the  normal.  The 
result  is  a  fourth  order  polynominal. 


a4Q4  +  OfjQ3  +  a2Q2  +  ai^  +  a0  =  o,  (9) 

where  aQ  -  (xg2  +  zg2 )  -  xg2  -  Zg2, 

“l  =  -2  <X0XS  +  zozs  "  xs2  "  zs2)* 

« 2  =  -(x  2  +  x  2  +  z  2  +  z  2  -  4x  x„  -  4z  z„  +  2  (x  2  +  z  2 )  (x  2  + 

*  o  s  o  s  OS  OS  o  O  S 

zs2>)» 

a3  =  “2  (xoxs  +  zozs  ■  xo2  -  zo2>» 

“4  ■  <xo2  +  zo2)  '  xo2  ■  zo2- 

The  equation  may  be  solved  by  synthetic  division  and  the  Newton-gradi¬ 
ent  method  or  other  numerical  solution  techniques.  (Hildebrand,  p.  451). 


There  are,  in  general,  four  solutions  to  equation  (9),  each 
corresponding  to  a  point  on  the  cylindrical  surface.  The  point  is  a 
reflected  image  of  the  light  source  if  the  light  source  and  pilot's 
position  are  on  the  same  side  of  the  surface.  This  is  the  case  if  the 
ratio  of  the  cosines  of  the  angles  between  the  rays  PgPr  and  PQPr  and  the 
surface  normal, 


Qso-  1  ~  anxs  ~  cnzs  (10) 
3  “  ~  cnzo 

is  larger  than  or  equal  to  zero;  i.e. ,  Qgo^0* 

The  computary  scheme  worked  out  above  for  the  simple  cylinder  case  is 
used  as  a  solution  space.  The  coordinates  of  the  positions  of  the  light 
source  and  pilot  are  transformed  into  their  equivalent  values  in  the 
solution  space.  The  reflection  points  are  computed  and  the  coordinates  are 
transformed  back  into  the  original  problem  space.  The  coordinate 
transformation  from  the  problem  space  to  the  solution  space  consists  of 
four  steps:  (Newman,  Ch.  3) 

1.  Translate  the  origin  (xc,yc,zc)  of  the  cylindrical  axis  into 
the  origin  of  the  coordinate  system  of  the  solution  space,  (2)  rotate  the 
cylindrical  axis  (ac,bc,<^)  clockwise  about  the  x  axis  until  it  is 
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aligned  with  the  y  axis  of  the  solution  space,  (3)  rotate  the  cylindrical 
axis  clockwise  about  the  z  axis  until  the  z  axes  are  aligned,  and  (4) 
scale  the  original  space  by  the  cylinder  radius,  R£. 

The  translation  of  the  coordinates  (x,y,z)  into  their  equivalent  value 
(xl,yl,z*)  in  the  solution  space  is  performed  by  matrix  multiplication; 
i.e. , 

(xl,y»,z»,l)  -  (x.y.z.l)  M-pMxMgMg.  The  corresponding  matrixes  are 
given  by: 

mt  -  /  1  0  0  0  \ 

(  0  1  0  0  \ 

l  0  0  1  0 

'  -*c  -yc  _zc  1  '  i 


0 

1 

0 

0 


0 

0 

1 

0 


where  v. 


4 


*c2 


+  Cc* 


Once  the  reflection  points  are  solved,  the  position  coordinates  (xrl  , 
y r * ,  zr*)  are  transformed  back  into  their  equivalent  values  (xr,  yr  ,  zr) 
in  the  problem  space«  This  operation  is  done  by  an  inverse  matrix 
transformation;  i.e.,  (xr,  yr,  zr,  1)  ■  (xrl,  yr* ,  zr!,  1)  M,“* 


M, 


M,-‘ 


M. 


Mf"1 .  The  corresponding  matrixes  are  given  by: 


1 


0 


0 


0 


1  0 

0  1 

yc  zc 


0 

0 

1 


» 


(12) 


0  0 
1  0 
0  1 
0  0 


0 


0 

0 

i/rc 


The  reflection  points  are  checked  against  the  window  edges  to  ensure 
enclosure  (Smyth,  July  1977,  App  A).  The  light  reflectance  is  computed 
from  the  angle  of  incidence  between  the  ray  and  the  surface  normal  (Smyth, 
June  1977,  App  A). 
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APPENDIX  C 

COMPUTER  PROGRAM 


Computer  Program 


The  computer  program  used  to  compute  the  transmittance  and  coordinate 
values  of  the  matching  sample  Is  attached  below.  The  listing  Includes  the 
program  and  subroutines,  the  raockup  frame  vertices  Input  data,  the 
position  and  orientation  of  the  mockup  on  the  test-area  floor,  and  the 
positions  of  the  light  source  used  In  the  test.  The  program  Is  written  In 
the  Fortran  IV  language  and  w*s  run  on  the  CDC  7500  Computer  System.  The 
subroutines  are  listed  as  follows: 

1.  CONTL-  Main  program  calls  for  read  In  of  data  and  computations. 


2.  READV-  Subroutine  reads  in  frame  vertice  coordinate  data. 
Called  by  CONTL. 

3.  NORML-  Subroutine  computes  surface  normal  for  each  canopy 
window.  Called  by  CONTL. 

4.  FLEG-  Subroutine  reads  in  mockup  position  and  orientation  on 
test-area  floor.  Called  by  CONTL. 

5.  SLIGP-  Subroutine  reads  in  test  light  position  on  test-area 
floor  and  calls  for  computation  of  light  images  for  the  left  and  right  eye 
of  the  observer.  Called  by  CONTL. 

6.  STPS-  Subroutine  converts  the  light  source  position  into  the 
coordinate  system  of  the  mockup,  and  calls  for  the  computation  of  the 
direct  view  and  primary  reflected  Images.  Called  by  SLIGP. 

7.  PRNT-  Subroutine  prints  out  direct  view  and  reflected  image 
point  information.  Called  by  SLIGP. 

8.  MTRNL- 
MROTX- 
MROTY- 
MR0TZ- 

MSCAL-  Subroutines  called  by  STPS  to  convert  light  source 
position  on  test-area  floor  into  coordinate  system  of  mockup. 

9.  C0MP0-  Subroutine  controls  computation  of  direct  view  image 
position  coordinates  and  transmittance.  Called  by  STPS. 

10.  C0MPL-  Subroutine  controls  computation  of  primary  reflection 
image  position  coordinates  and  transmittance.  Called  by  STPS. 

11.  INTEC-  Subroutine  computes  direct  view  image  point  on  a  planar 
surface  for  a  given  light  source  position  and  checks  whether  point  is 
enclosed  by  surface  edges.  Called  by  COMPO  and  COMPL. 
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12.  INTCY-  Subroutine  computes  direct  view  image  point  on 
cylindrical  surface  and  determines  whether  point  is  enclosed  by  edges  of 
cylindrical  window.  Called  by  COMPO  and  COMPL. 

13.  TRSCY-  Subroutine  converts  coordinates  of  a  cylindrical 
coordinate  system  into  a  rectangular  system.  Called  by  INTCY  and  INTS  when 
testing  a  surface  point  against  surface  edges. 

1A.  INTPL-  Subroutine  computes  primary  reflection  point  coordinates 
for  a  planar  surface  and  determines  whether  the  point  is  within  the 
surface  edges.  Called  by  COMPO  and  COMPL. 

15.  INTCL-  Subroutine  computes  primary  reflection  point  coordinates 
for  a  cylindrical  surface  and  determines  whether  the  point  is  within  the 
surface  edges.  Called  by  COMPO  and  COMPL. 

16.  TRANC-  Subroutine  converts  point  in.:o  cylindrical  solution 
space.  Called  by  INTCL. 

17.  SPOLY-  Subroutine  computes  coordinates  of  image  points  in 
cylindrical  solution  space.  Called  by  INTCL. 

16.  SOLN-  Subroutine  checks  image  point  to  ensure  that  it  is 
reflected  and  controls  conversion  to  problem  space.  Called  by  INTCL. 

19.  INUC-  Subroutine  converts  image  point  from  cylindrical  solution 
space  into  problem  space.  Called  by  SOLN. 

20.  INTS-  Subroutine  determines  whether  reflected  image  point  on 
cylindrical  surface  is  enclosed  within  surface  edge.  Called  by  INTCL. 

21.  COMP-  Subroutine  computes  transmittance  and  reflectance  for 
image  po^nt.  Called  by  COMPO  and  COMPL. 
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SNYTF>. S T H |  l , 

ACC  OUNT  »  ME •*  *  . 
f IN  t  SL  #  P  *  . 

MAP, ON, 

LGO. 

PROGRAM  CONTI  ( INPUT  Pul  .  1  AP  £2-  INPUT.  TAPE  5-1)1  T  P  0  T  ) 

CALL  REAOV 
CALL  NO»ML 
CALL  ELEG 
CALL  SHOP 
STOP 
FNC 

SUBRLUT  INI  REAOV 
C  READ  IN  iUREACF.  VERT  i C i  S 

COMMON 7C AN/NT.NB*  NC  .NP.NA.NO.NV ( 1 00 ) * P X V ( I  00#  c ) . P Y V (  1 00 ,  t  I  , P Z V 1  1  OC 
0*6  ) 

CUMPl.N/  Vt  B  I  /  XV  {  2 CO)  »n  1  7  00)  .  ZV  (  2  00)  »NVR  (  100*8  ) 

COMMrNACYL/NCY.NSCl 1C  )  *  N  S  P  ( 10* 10)*XC(10).YC(10>»ZC( 1  C )  ,  A) 1  10) *  C  E  1 1 
OCI »  CF ( 10 ) *  RC ( 10) 

C  C  M  f-  C N/PIim  /  Xii-VO*  ZL* 

RtA{;(7»1000) 

iccc  Frtutin 

R  F.  A  r  (  ? »  1000) 

R  E  A  D  (  2  *  1001  )NT*NE*NC»NE*NA*NO 
1001  EORMA  f  <61  it,  I  3  )  ) 

REATl’.lOOO) 

RE  AT (2,  IOC?  )(NVti)*I»l*N,)| 

10CL  FORMAT < A< 2X*  I  3  )  ) 

B  c  A  C  (  2  *  1000) 

nn  ic  j-i.n 

REA0(2,1002)(NVP(1.JI.1«1.N0) 

10  CONTINUE 

REA012. 1000) 

BEAUS.  1003)  1  AV(  i  )*  I-l.NT) 

1C03  E0RNAT(a(7X,) 7.3) ) 

READ!?, 1C00) 

RE  ADI  2,  1003  )  ( Y V ( I ) .  I-  1. NT  ) 

READ! 7,1000) 

REACM2. 1003  )  1  ZV<  1  )  .  I-1»NT  ) 

REA 0(2*10 00) 

RE  AO ( 7*  1002JNCY 

IE  INC  Y  .EO  .0  JGUTC  2b 

RE  A(.  (2.  1002)  (N!t  (  1  )  » I-l.NCY  ) 

00  20  I  ■  1  »  NC  Y 
KP-NSC  1  1  ) 

20  READ(  2.1002  KNbPU*'<)*K-l*KP) 

READ(2»1004)(XC(I)*YC(I).ZC(II.AF(l).cE(l).CE(H»PCn)»l-l.NCY) 
10GC  F0RMAT(7(7X»E7.3) ) 

2b  CONTINUE 

PE  Af  1  ■>,  1000) 

BEAC(2*  1003  Mn.  Yl»ZO 

RETURN 

END 

SUBROUTINE  NOR  M  L 

C  establile  surface  mirral  f op  each  plate  surface 

C  SURFACE  NJP"ALS  DIRECTED  TOrARD  COCKPIT  INTERIOR 

COPHON/C  AN/NT ,NB .NC* NP.NA.NO* NV (I 00).EXV(100.H).PYV(100.t)).PZV(100 
0.6  ) 

CUPNLN  /  VERT/XV(?00)»YV(700)fZV(2  00)»NVR(100*’F) 


mu  Tint  I  *  -  v  r  tty  ntiAMOA^iUM 

■m  cu?  ■ 


CUMMGN/NQRM/aXN(100)»AYN!1QO)>A7N(1001 

CJMMCN/CYL/NCY,NSC!lO)#NSP<lO»lO>,Xu<lC;),YC(10)»ZC<lO),AEClO)»BE(l 
QOlfCi  ( 1,0  )  ,  RC  (  10  ) 

COMMON/ P 11 CT / XO, YO# ZO 
DU  1C  1  •  1  #  Nf5 
NK»NV< I  ) 

Dll  i  K  •  1 ,  NK 
KV«NVR ( I »K ) 

PXV  < I »K  )  aXV!K V  ) 

PYVII.K)-YV(KV) 

PZV(I»K)»ZV(KV) 

5  CONTINUE 
K  ■  2 

7  A 1 • P  X  V ( I »  K  )  -P  X  V  (  I.  1) 
az«pxv(  i »«♦ n-pxv  ( i  .1 ) 

Bl-PYV! I»K)-PYV( I  .1  ) 

B2«PYVU»K41)-PYV(I»1I 

ci»pzvi  i»k  i-pzvn»i  > 
cz-p/vi  i»K»n-pzvu , i ) 

PlaS<jRT(Al**i>Bl**2*Cl**2) 

P2»SCkT(A?**0+B2*«24C2**2) 

A-(Ai*A2+bl*B?+Cl*C2)/<Pl*P2> 

1F(AES( Al.LT.l.  I  W  TO  ? 

K  «  K  ♦  1 

IF (k.EU.NR  )  GO  TO  10 
GO  TO  7 
9  AN  ■  AC  OS  (  A  ) 

R«l./<Pl*P2*SIMAM  ) 

AXN(I)a-(Bl*Ci~Cl*P2I*R 
A YN ( 1 )■♦<  A1*C2-L1+A2)*R 
A  2  N ( 1 J  «-( A1+D2-A2+B1  >*R 
10  COM  I  SUE 
RETURN 
END 

SUBROUTINE  FLFG 
C  CRAFT  POSITION  IN  SCtNE 

C  CRAFT  ORIENTATION.  DIRECTIONAL  COSINES  OP  Y-AXIS  AND  Z-AXIS  OF  CRAFT  l 
C  SCtNE  SPACE.  T- AX  IS  CIRECTEO  BACKWAkTS  ALONG  LONGITUDINAL  AXIS.  Z-Axl., 
C  TOWARD  VERTICAL 

COMMON/ SC /XL. YL./L, Al.ttL.CL, AAL.BBl.CCL 
C  CRAFT  POSITION  AND  ORIENTATION  IN  SCENE 
Pt AC(2.999) 

999  FORMAT!) 

REAU2,100A)XL,U,ZL.AL.BL.CL.AAL.BBL,CCL 
10G  A  FORM  A  T(2X#i!FK.<.»2X)/2X,3(Fl0.^»2XI/2X.B(F10. 4.2X1) 

WRITE (3.1 OC  5 ) 

100b  F  OR  MA  K2X,  •  POSITION’  /AX,  •*-POS,»<»X»»Y-POS»»<tX.,Z-POS'l 
WRIT E< 3,1 006)XL»YL»Zl»AL.BL,CL»AAL»BBL»CCL 
1006  F(1PMAT(AX,3(F10.<I,2X)/<,X,3(F7.A,2X)/^X,3(F7.4,2X)) 

RETORN 
END 

SUBRCUTINL  SLIGP 
COMNr  N/GT/X, Y.  Z 
COMMCN/P ILJT/XO, Y J, ZC 
DATA  DE/2./ 

DATA  AP/1H-/ 

XOO»XO 

RF  A U (  2»  1  GC  1 ) ND  P 
WR I TE ( 3» 1 001  »NCP 


t. 


"t 
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1001  FORM*!  l?x,13) 

P0  10  I  • 1  *  NOP 
i>tAP(?,1002)X»Y,2 

1002  FCRMAn*l2X,F10.A)) 

WRITE(3»IC00)1,X»Y#Z 

1000  FORMTItH  ,  •LC.«»»lS»2X»'Xl-*,ri0.A»2X,'YL«',f  10  .  A  ,  2 X , • Z L • • » f 10 . A  ) 
XO-XOO-3. 

00  7  J  » 1 »  2 
XO-XO+OE 

WR  1TEI  3. 10031  X0, YO»ZO 

1003  FORMAT <2X,«filuT-LY£  F  OS  [  J  1  UN»  ,  3  (  2  X,  F  1  C  .  A  1  1 
CALL  STPS 

CALI  RRNT 
7  CONTINUE 

WRITE  (  3,  l  004  MAP,  K-  1,80) 

ICC*  FPRMT(2X,flOAl) 

10  CONTINUE 
RETURN 
END 

SUBROUTINE  STPS 

C  COMPUTES  AND  STORES  STATIC  ANALYSIS  CF  LIGHT  SOURCE  RtFLtCT  1LNS 
COMMUN/LIUHT/XL.YL.ZL 

C0NECN/SC/XS,YS,ZS,AS,6S,CS,AAS,88S,CCS 

COMNUN/GT/XG, yg,zg 

COMKiN/FP/NSL»XP,Y?,ZP,TP 

X«XS 

Y-YS 

Z«ZS 

A1«AS 

Bl-bS 

Cl-CS 

A  A 1  ■  A  A  S 

88  l«tJ8S 

CC  1  *CC  3 

VC«SCRT(Bl**2+Cl+*2> 

CALL  M  8  G  T  X  (  81/VC, C 1/VC, AA1, 381, CC1  ) 

CALL  MP0T2(vC#Al,AAl,oPl»C:i) 

XL«XG 

YL-YG 

Zl«ZG 

CALL  MTRNL<-X,-Y,-Z,XL, YL, ZL) 

CALL  NR0TX(61/VC,C1/VC»XL,YL,ZL 1 
CALL  MR0TZ(VC,A1,XL,YL,ZL) 

CALL  NR.3TY«CCl#AAl,Xl#YL,ZU 
CALL  COMPO(NSL,XP»YP,ZP, TP) 

CALL  CONPL 

RETURN 

END 

SUBROUTINE  PRNT 

C  PRINTS  01 T  LIGHT  SOURCE  DATA,  DIRECT  VIEW  AND  REFLECTION  POINTS 
CQMHr.'N/LIGHT/XL,YL»ZL 
COMMON/ FP/NSR,XP, YP,ZP»  TP 

C0MMCN/FY/NSl,li(10),XS(l0),YStl0),ZSll0l,lR(L0),XR(10)»YR(lC),ZR( 
010 ) , TR ( 10  ) 

WRITE)3,1000)XL,YL,ZL 

loco  FGRRATUH  ,ZX, 'LIGHT  SOURCE • / 2X , • XL • • » U  0.  A , 2 X.  • Yl • « , U 0. A, ? A,  •  Z L » 
0 • , F 10«  A ) 

IF (NSP.EO.O)GQTO  5 
WR1TE(3,1001)NSP,XP,YR,ZP,1P 


mTS  flSt  IS  8*8? 

"  '  •  *  •  •*>  \  . 


WWtTT  r.TASt-T 


AO 


1001  FORMAT (IH  ,2X,*0IPtCT  V I tM • /2X. I 3, 2X, 4 ( 3X. F 10 . 4 t) 

!>  CONTINUE 

IF<NSL.E,0.0)RETUPN 
OG  E  l  ■  1 ,  N  S  l 

fc  WRITE  <  3»1C02>  IS(I),  XS(  I».  YS(  I  ),ZS<  I ),  IR(  I).  XR  (  I  ),  YR<  DiZHIlli  TP  (  I) 

1002  FOPMATdH  ,2X, 'REFLECTION  P(J  1NT  •  /  (  2  X  ,  2  (  1  2  »  2  X  ,  3  (  F  10 . 4  »  2X  )  )  »  F 1G  .  4  )  > 
RETURM 

END 

SUBKf UT  INF  CC.N.POINSL,  XP#YP»ZP>  TP  ) 

C  COMPUTES  IMTLFCrPlIOK  POINT  OF  STRAIGHT  LINE  FAY  BtUtEN  PILOT  AND  SOL 
CCMM.CN/CAN/NT,NB,NC»NP»f:A,ND.NV{lCC)  ,PXV(  100,  8  I  »T  TV(  100,8  ),PZVUGt 
0,8) 

COMMON /LI GHT/XL, YL, ZL 
COMMCN/PILOT/XO* YO.ZO 

CGMMDN/LINE/AS»HS»CS»XS»YS»ZS»tC»BC»CC 

NS  L  *C 

TPM. 

XS«XC 
YS«YC 
Z  S  *ZC 

R»SOFT(<XL-XO)**2+<YL-YO)**2+<ZL-ZO)**2) 

ASMXL-XOI/R 
BS» (YL-YC  )/R 
CS«(ZL-ZO)/R 
DO  1C  I S *  1 »  NO 
ISK  *0 

IF  <  IS.LE.NPKALl  INTECUSK,  IS,X?,YP,ZP> 

1F(  1S.E0.NA)  CALL  INTEC(ISK,1S,XP,YP,ZP) 

IF  (  JS.GT.NPK  ALL  I  N  TC  Y  ( I  SK  ,  I  S  ,  XP  »  Y  P  ,  ZP  ) 

IF  (  ISk.GT  •  u  )  G  D  T  f.  20 
1C  CONTINUE 
RETURN 
20  CONTINUE 

IE ( JS.LE.NC  jKfcTU-N 
IF  I IS.EO.NAIGOTU  30 
NSL*  ISK 

CALI  COMP ( ANG.KT, TT  ) 

TP-TFMT 

RETURN 
30  CONTINUE 

CALL  COMP  I AN&, R  T,  TT  ) 

tp-tfmt 

GO  TP  10 
END 

SUBROUT  INE  CUMPL 

C  COMPUTE  PRIMARY  REFLECTION  POINTS  FOR  SINGLE  SOURCL  OF  LIGHT 

COMMON /CAN/NT, NB,NC,NP,NA, NO, NV(100),PXV1100»E),PYV(1GO, 8), PZV1100 
0,8) 

COMMUN/NURM/AXN ( 100) , AYN( 100), AZN( 100) 

C0MM0N/CYl/NCY,NSCI10),NSPC10, 10 ) > XC ( 10 ) , YC < 10 > , ZC ( 10 ) » AE ( 10 ) , B E ( 1 
00) ,CE  < 1 0) ,«C( 10) 

COM.MCN/LINE/AS,BS,CS,XS,  YS,  ZS,  AC,  BC,CC 
COMMON/LIGHT/XL, VL,Zl 

COMMGN/CPFF /NN, AN ( 10 ) , SS ( 10 ) , NS, XRR 1 10 ) , YRR ( 10 ) »  ZRR { 1C ) 
C0MM[N/FY/NSL»ISX<10)»XSX<10I,VSXC  101,  ZSX<  10),  1<?X<  101,  XPX<  1C  ),»*»< 
010),ZRX(10),TRX(1C) 

NSL-C 

NSS-MC+1 

OC  2C  IS"NSS,ND 


IF  US.LE.NP1G0I0  30 
1 F  (  1S.GT.NPIGCTU  NO 
60  TC  20 
3G  CONTINUE 
AC«AXN< IS ) 

BC«AYN( IS  » 

CC-A7NUS  I 

*s»p*vns»n 

ys»pyvus.ii 

ZS*P7V(1S#1) 

ISK»C 

CALI  INTPLUSK.  1 !.  I 
1M  ISK.tU.OIGOTO  20 
CALL  C  I'MP  ( ANG  »  R  *  T  T  ) 

00  32  I ■ 1 » Nn 
ISK-C 

IM  I.tO.ISIGOTO  32 

IF  U.Lfc  .NF1CALL  INTEC  C ISK, l»XR»YR#Z&  > 

IM1.GT.NPICAIL  lMCY(ISK,I»XFi  YK,2») 

IF  ( ISK.LE •OIGUTC  32 
IF  (I.LE  .NCJGOTO  20 
IF ( I.EO.F AJGUTU  36 
RX»AS*AC«bS*0C  +  C.S*CC 
AS«AS-2.*RX*AC 
BS«bS-?.*PX*0C 
CS*C  S-2  « *fc  X *CC 
CALL  COMP (ANG »KT»TI  ) 

T  ■  P  ♦  1  T 
NSL«NSL*1 

WRITE(3ilC01>NSt»I5.>5»YS#ZS#l»J<F#Yt«»Z*»T 
1001  F0*KAT(2X#IS»2X#I3»?X>3(F10.N»2X|#2>#13#2X,N<F10.4»2X>) 
ISX(NjL)*IS 
XSXCNSL  »»XS 
YSX(KjL)*YS 
ZSX(NSL)«ZS 
IRXIKSL3-1 
XP  XiNSL  )*XK 
YRX (NSL ) • YR 
ZRX (NSL  1  «ZP 
TP  x (Nit  )■! 

GO  TO  2C 
3fc  COM iNUt 

CALL  CUMP(ANG»RT.TT) 

R*MTT 
32  CONTlNUfc 
GO  TO  20 
NO  COM  IN  IE 

00  NU  I C • 1 »  NC  Y 
KP-NSC (  IC  I 
DO  Nil  KC*1#KP 
I f  C  NS  P ( IC  »KC) . c  0  .  i S  )  GO  TC  N2 
GO  TT  Nil 
N2  CONTINUE 

CALL  INTCLUC1 


IF (NS.GT.CICALL  INTS(IC.IS) 


IF (NS.EO.OIGUTO  20 
00  NN  IK* 1*  NS 
Xi «XpK ( IK > 

YS«YRR(  IK  I 


m  t’CT  m ^ .V*"  ^ 

rot*  cwi  »-»*v 
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7S«ZkR(IK  » 

RL*SORT((XS-XL  YS-YU**2MZS-ZL  J**2  I 

AS*  (Xl-XSWRL 
as- tYL-YS  »/PL 
CS»(ZL-ZSI/RL 

KO-IXS-XC  < I C ) )*AE ( I C I MYS-YCL I C » I *BE ( IC I ♦ ( Z S-ZC ( IC I  )  *CE ( IC  J 
AC«<XC(ICI-XS*KO*AC(IC»»/ftC(IC) 

BC» CYCLIC  >-YS*kO*BFCIC»  J/RC(IC> 

CC*CZC<IC)-ZS«-RC*CE(ICI)/RCCIC) 

CALL  C  LMP ( ANG  »  k  »  T  T  ) 

DO  A3  I  “  1 »  NO 
ISK*C 

I F ( 1 .£  Q. I S ) GOTO  A3 

IF ( 1  .LE «NP  )CALL  I  NT  EC  < I $K . 1 , XR  ,  Yk  ,  ZR  ) 

IFU.GT.NPICALL  J  NTCY  { I  SK,  I »  XR  ,  YR  ,  Zk  ) 

IF ( ISK.LE .OIGOTO  A3 
A 5  CONTINUE 

IF (  1  . LE . NC ) GOTO  AA 
I F ( I  .tO.NAIGUTC  Aft 
PX*A S*AC+PS*8C+CS*CC 
AS«AS-2.*RX*AC 
GS*PS-2.*RX*BC 
CS«CS-2 .♦Rx*CC 
CALL  C QMP ( ANG  »  R  T»  T  T  ) 

T»k*TT 

NSL*NSL*1 

WRIT(<3.1001)NSL.IS.XS,YS»Zj»l.XP»|R»ZK>l 
GO  TO  AA 
A6  CONTINUE 

CALL  COMP ( ANG.RT.TT  ) 

R«R*TT 
A3  COMINUE 
AA  CONTINUE 
GO  TO  20 

All  continue 

20  CONTINUE 
RETURN 
ENC. 

SUBROUTINE  INTPLI  I SK. I S ) 

C  COMPUTES  PRIMARY  REFLECTION  POINT  OF  LIGHT  GIVEN  POSITION  WITH  PLANE  G 
C  VERTEX  ANO  SURFACE  NORMAL  ANO  POSITION  OF  VIEWER 
C  DETERMINES  IF  KAY  STRIKES  CONVEX  SURFACE 

CCMPCN/CAN/N' T»N6»NC»NP»NA»ND»NV(100).PXV(10Q»8)»PYV(1C0.8).PZV(100 

0.8  I 

COMRON/UNE /AS.  BS.CS.XS.YS.ZS.AC.BC.CC 
COMMCN/LIGHT/XL.YL.ZL 
COMMCN/PILOT/XO.YO.ZO 
QL«/C*(XL-XS>*BC+(YL-YS)ACCMiL-2S) 

00«AC*(X0-XS)+8C*< YO-YS) tCC-M  ZO-ZS  J 
ALF-OL/OC 

IFIALF.LT. O.IRtTURN 

C  LIGHT  SOURCE  ANO  PILOT  ON  SAME  SICE  OF  PANEL 

RLO*SQRT< (XL-X0»**2*( YL-YO I **2 ♦ l ZL -ZO I  * *2 » 

AL0«(XO-XL)  /RLO 
BLO«(YO-YL)/RLO 
CL  0* ( ZO-Z  L I /R  LC 
AV«ALF*RLC/ C1.*ALF  » 

XE»XLAALO*AV 
YE ■YL+8L0*AV 


2  E ■ 2 1 ♦CLP* AV 

RE*AC*CXE-XS>+BL*tYF-YS)»CC*(2l-ZS) 

Xi«XE-AC»*R£ 

YS«Yf-BC*R£ 

ZS«ZE-CC*R6 

Rl»S«RT<  tXS-XLJ**2MYS-YU«'*2P<ZS-Zi  »**2> 

AS»(Xt-XSI/RL 
BS-tYL-YS I /fit 
cs*m-zs  I /RL 
I N  «KV (  IS) 

CO  10  I  •  l  »  1  N 

IC - !♦  I 

IF  ( I.EQ.  1M  IC«1 
A 1  ■  P  X  VC I S  » I >-X0 
81 «P V V  C IS*I)-YO 
Cl-PxV( IS/1 )-Z0 
A2-PXVC  Ij.IO-XO 
B2-PYVC IS. IC)-YC 
CZ‘FZV<lSt  10-20 
Pl-XS-XO 
P2-YS-Y0 
P3»Z$-Z0 

0-Pl*(3l*C2-B2*Cl)-t»?*<Al*C2»Cl*A2>  +  P3A(Al*82-Bl*A2) 

IF(G. IT. C. RETURN 

C  RAY  STRIKES  SURFACE  UN  ENCLOSED  S10E  OF  SURFACE  EDGE 
10  CONTINUE 

C  RAY  STRIXtS  INClOStC  SURFACE 
ISK-1 
RETURN 
END 

SUBROUTINE  INTCUId 

C  OETERMlNt  CYLINDRICAL  REFLtCTI )N  POINTS 
CCNPLN/CICHT/XLL, YLL,ZLL 
COrMCN/PILOT/XOO, YCO»ZOO 

CO(*Nf;N/COFF/NN,  AM  1  0  I  ,  CS  <  1C  >  »  K  S*  XRR  (10 >  *  YRR  110 )  ,  ZRR  I  10  J 

CONhCN/Pl AC/XL, YL  »ZL»XO» YO«ZO 

DATA  NN  /  W 

XL ■ XL  l 

YL-YLL 

2L-2LL 

CALL  TRANCUC»XL*YL»7l  ) 

XO-XOO 

YO«YOO 

ZO-ZOO 

CALL  TRANCI IC#XO, YO.ZO) 

ROL«!ORn<XL-XO»**2*<YL-YO)**2HZl-ZO>**?> 

BOLMYO-YLIZROL 

AN  Cl >•< X0**2  +  Z0P*?>**2-X0**2-Z0**2 
AN  (?) •-£.*( XO*XL* 20*21 -X0**2-Z0**2> 

ANC3)»-<XC**2nXL**2+ZO**2nZL**2-A.*XO*XL-4.*ZO*ZL*2.*(XOA*?*ZO**2) 
Q*«  Xl**2*7L**2  I  ) 

ANCA)«-2.*( X0*XL*70*ZL-XL**2-ZL**2I 
AN(5I«CXL**2+ZL*A2»**2-XL**2-ZL**2 
CALL  SPOLY 
N$»0 

DO  10  1 • 1 #NN 
AF-RCL/CCS1 I)*l.  ) 

RE-SORTCC  XO  +  CSCI )*Xl )•♦?♦! Z 0»CS C I) *Z L » ♦ *2 1 
AtMXO»CSm*Xt)/KE 


ce-izd+csumzu/re 

YE  ■  Yl  ♦  AM  60L 
XR«Ak 
YR-YE  ’ 

ZR  *CE 

CALI  SOLN ( I  C • XR • YR  »  ZR  I 

XP«-AE 

YR-YE 

ZR»-CE 

CAL l  SDLN(1C.XR»YR#ZP) 

10  CONTINUE 
RETURN 
END 

SUBPCUT INE  SOLN(  IC» XR.YR,ZRI 

CQMMlN/CCFF/NN>AM10)»CS<10»»NS»XRM10l»YP«{10>»ZRR<iC> 

COKKLN/PLACm#YL»ZL#XO»YO,ZO 

KL ■ SOR  T ( (Xfi-Xl )**Z+IYR-YL )**2MZF-iLl**2) 

RC-SCRTI  I  XR-X  £)>♦♦?♦  ( YR-YO I *♦ 2* ( ZP-Z 0 > ** ?  | 

ol-( i.-(xs*xi-*zr*zi  n  /pi 
00«O.-(XR*X0  +  ZR*Z0))/P0 
IFICL/OO.LT.O. IRETURN 
Q-Cl+OO 

IF (C  .EO.O . J  RETURN 

IF  «  ft  e  >  t  <  OL-QO » / 0  I  .OT.. 01 (RETURN 

NS -N  S ♦ 1 

CALL  INVC (  ICtXR, YR, Z«) 

XRP ( N3 ) ■ X  R 
YRP  INS  I -YR 
ZRRINSI-ZR 
RETLPN 
END 

SUBRC'JTINE  MTRNUXT,  YT/ZT#X,  Y,Z> 

C  TRANSLATION  BY  POSITIVE  SCALARS 
X-X+XT 
Y-Y+YT 
Z-Z+ZT 
RETURN 

end 

SUBROUTINE  NROTX(A,B»X»Y»Z) 

C  ROTATION  ABOUT  POSITIVE  X-AXIS  OF  Y-AXIS  IN  POSITIVE  DIRECTION 

A  X  «  A 

BX— E 

YV«AX*Y-BX*Z 

ZV-BXAY+AXAZ 

Y«Y  V 

Z-ZV 

RETURN 

END 

SUBPUUT INE  MROTYIAt 8»X,Y#ZI 

C  POTATION  ABOUT  POSITIVE  Y-AXIS  OF  Z-AXlS  IN  CLOCKWISE  DIRECTION 
AY-A 
BY  •-P 

XV«AY*XfBY*Z 

ZV«-BY*X4AY*Z 

X*XV 

Z-ZV 

RETURN 

ENO 

SUBROUTINE  MROT Z (A , B , X , Y » Z  I 


C  POTATION  ABOUT  POSUlVt  2-AXIS  uf  Y-AXIS  IN  CLOCKWISE  DIRECTION 

AZ  *A 

BZ  »  P 

XV«AZ*X-BZ*Y 
YV»BZ*X*AZ*Y 
X*  X  V 

y«yv 

RETURN 

END 

SUBROUTINE  MSCAL ( A,*»Y.Z  ) 

C  M AGN IF  I C  A  T  l  ON  8Y  PLXlUVf  SCALARS 
X«  A** 

Y»  A*  Y 
Z«A*2 
RETURN 
END 

SUBRtlUT  IN  t  TRANC(IC«X,Y,Z) 

C  TRANSFER  INTO  CYUNCR1CAL  SuLUTION  SPACE 

C0NMCN/CYL/iNCY,NSC(10)#NSP(10*10l»XC(10I»YC(l0I»ZC(I0l*AC(10)fBC(I 
OOI »CC (  101 »»C< 101 

CALL  MTRNL(-XCIIC*»-VCnCl»-ZC(IC)#K#Y»ZI 
VC«SCAT(tJCUC)**2*CClICJ**2» 

CALL  MR0TX(3C<IC)/VC>CC(ICWVC»X»Y»Z» 

CALL  MROTZTVC.ACTIC  I.X,Y,Z> 

CALL  MSCAKl./KCUCIf  X.Y.Z) 

KETUFN 

END 

SUBFLUTINE  INVCUC.  «•  Yi  L  ) 

C  INVERSE  TRANSFORM  FRC.M  CYL INORICAL  SOLUTION  S*>ACE 

COMMLN/CU/NCY»NSCnO>»NSPUO»10)#XC(10>»YCUO>»ZCUO>*ACUO>»8C(l 
00)ft(  UO).RC(IO) 

CALL  NSCAL1RCUC)»X»Y,Z) 

VC-SCRTTBCTIC  »**t'*CC(  IC>**2> 

CALI  MPijTZI  VC*-ACnCI»X,Y#Z) 

CALL  MRQTX(BC(ICI/VC.-CC< ici/vc#a# Y,Z) 

CALL  NTRNLIXC  I1C)»YCUCI  »  ZC  UC  I  »X»  Y»  Z  > 

RETURN 

END 

SUBRUU  TINE  SPQLY 

C  REAL  PLOTS  TO  POLYNOMIAL  EQUATION  USING  SYNTHtTIC  DIVISION  AND  NEWTON* 
C  GRADIENT  METHOD.  ASSUMPTION  THAT  NO  COMPLEX  ROOTS  PRfStNT. 

CONMLN/CUFF/NN, AN  ( 10  >  »  XS  <  10  >  >  N  S»  XRR  <  10  >  *  YRM1  0  >  »  ZR  R  I  10) 

CIMFNSIuN  BN ( 10) 

N«NN 

00  3L  I • 1 #  NN 
1MAM1I.E  J.O.IGUTO  41 
1MaMZ>«E0.0.  >60 TO  42 
X£«-ANI1> / AN  <  2  * 

GC  TC  5 

41  XF«G. 

GO  TC  5 

42  IF ( AN ( 3  I  •  EW  •  0 • ) GOTO  43 
»c**CRT(ABS(AN(I)/AN(3>>> 

CO  TO  5 

43  IF  IAM4).E3.0.»G0TU  44 
XE-(-ANU  J/ANI4  >  >**.33333 
GO  TC  5 

44  XE*(ABS(AN(1>/AN(S>)>*«.2S 
!>  CONTINUE 
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DO  20  J«l,100 
RNIN+l  l-AMN  +  1  ) 

C«0. 

00  10  K-’l.N 
L»N»1-K 

SN(L)"AN(L)*ltE*BN(L*l> 

C»C*XE+QNCL*1 ) 

1C  CONTINUE 

IF (C.cO.O.)GUTO  40 
/C 

CO  TC  15 
40  CONTINUE 
XR--10000. 

15  IF ( t  X  R  ♦  X  E l.EQ.O.IGOTO  20 

IF  CABS t «*»-*£)/ CXR*<E I). IT.. 00 II  GO TO  2  2 
XE  «XR 

20  CONTINUE 
22  XSCIWE 
N  «  N  —  1 

AN<N*1 )*8NIN*2) 

DO  3C  J  *  1  *  N 

L-N+l-J 

AN ( L ) *  BN ( L ♦ 1 ) 

30  CONTINUE 
PE  TCP N 
FNO 

SLBPLuTINE  INTSCIC.ISI 

C  DETERMINE  WHETHER  kfcHECTEO  RAY  POINT  IS  -ITH1N  ENCLOSED  SURFACE  UN 
C  CYLINDER 

CUPNON/C  AN/NT,  Nts,  NC  *  NP#  NA»  NO,  N  V  ( 100  )»  P  X  V  ( 1 00, 8  )  ,P  YV  (  100»fc),PZV(lC0 

0,8  ) 

CCMPON/COFF /NN, AN [ 1 0 ) , C S « 1 0 ) , NS, X S ( 1C  I » Y S ( 1 0  I , 2 S ( 10 > 

CUMKCN/PUnT/XO,YO,ZO 

XSS-XO 

YS  S ■ YO 

ZSS-ZO 

CALL  TRSCYCIC.XSS.YSS.ZSS) 

ISK-C 

00  2L  IK-1, NS 
XRP-XS  <  IK  1 
YRP-YSC IK  ) 

ZRR-ZS  C  IK  ) 

CALL  TRSCYCIC,XRR,YRR,ZRR) 

Pl-XRR-XSS 

P2-YRP-YSS 

P3«?PR-ZSS 

xi«pxvci£*n 

Yl-PYVC IS,1) 

Zl-PZVC  IS,1> 

CALL  TRSCYC  IC,X1,Y1  ,Z1 1 

A1-X1-XSS 

B1 ■ Yl-YSS 

C1-Z1-ZSS 

IN-NV ( I S ) 

DO  U  1 1  *  1 ,  I N 
I-IItl 

IFCII.EQ.INlI-1 

x2-pxvus,n 

Y2»FYvm,n  ;  ' 
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Z2-PZVI  IS,l) 

CALI  TRSC YI  IC» X2* Y2,Z2  ) 

A2«X?-XSS 

B2-Y2-YVS 

C2«Z?-ZSS 

0»Pl*<fil*C2-d2*Cl  >-l»2*<Al*C2-Cl*A2)  ♦P3*Ul**i2-Bl*A2» 

If IQ.LT.O. >GOTO  ?0 

C  PAY  STRIKES  SURFACE  ON  ENCLOSED  SIDE  OF  SURFACE  EDGE 
A1«AZ 
81»B2 
C1-C2 

10  CONTINUE 

C  RAY  STRIKES  ENCLOSED  SURFACE 

I$K»1SK*1 
XSUSKI-KSUK  » 

YSI  1SK)«Y$(  IK  ) 

ZS(ISK)«7$<  Ik) 

2C  CONTINUE 
NS  *  1 SK 
RETURN 
ENO 

SUBROUTINE  INTCYl iSK, ! S, X«, YR, ZR  ) 

C  CONFUTES  INTERSECTION  POINT  OF  LINE  WITH  CYLINDER 

CONNCN/CAN/NTtNo*  NC  *  NP»  NA,  ND  »N  V  1 1 OC  >  *  PX  V  UOO*  B  )  *P  Y  V I  1G0*  F  >  *  P  Z  V I  1  OC 

0,B) 

CQNNCN/LINE /AS,RS*CS,XS*YS,ZS»AC*6C»CC 

COFMCN/CYl/NCY,NSC( IOI.NSPI 10. 10)*  XC 1 10), YC I  10) »ZC ( 1C), AE (10 >,02 ( l 
Q0),Cl<l0>,HC<10) 

COrPC*N/PItOT/XP,YP,ZP 

C  DETERMINE  C YL INOEK ICAL  SURFACE  kHJCH  CONVEX  SURFACE  IS  A  PART  OF 
00  2  I«1,NCY 
KP-NSCU) 

CO  2  K •  I  *  K  A 

IF<NSP(I»K).Ea.rSI6nT0  5 
2  CONTINUE 
RETURN 
5  CUNTINUE 

C  DETERMINE  INTERSECTION  POINT 

xo«xcm 

YO- YC  1 1  ) 

ZO-ZCIi) 

RO-RCII  ) 

AC«AE  < I ) 

60«Pt ( I ) 

COCE  II) 

RCS«S<JRT<<XS-XO)**2HYS-YO>**2MZS-ZC>**2) 

ACO-IXS-XO/ROS 

BOC» I YS-YC ) /R  OS 

COC»(ZS-?C)/RCS 

A« AO*AstF C*HS+CO»CS 

B«RCS*( AO*AOO+BU*BOO+CO*COO> 

A1«AS-A*I (A0)A*2) 

A2*B  S-A* ( (bO)**2) 

A3«CS-A*( I CC) **£ I 
81»Al3*RuS-B*I  I  AC)**2) 

02«BOO*ROS-B*UBO)*A2) 

B3»C00*RUS“d*((CC)*A2) 

A»A1**2*A2*a2rA3*A2 

B«Ai*rtl»A2*o24A3»e3 
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C»Bl**?'»B2**2*b3R*2 

AB--E/A 

B8»  I  S  'JR  T,l  8*  *2-A  * (  C-R0**2  )  )  )  /A 
RS-AB-8B 

If  (PS.lT.C. IRS-AB+BB 
XR«X$*AS*KS 
ZK»Zb*Ci*RS 
YR«YS*BS*PS 

R«RS*( A0»AS»80*8S*C0*CS  » ♦ K 0 S  *  C AO*AOO  +  bO*BOO  +  CO* CCO ) 

Xl-XO*AC*R 
Yl- YC ♦GO^R 

zwoco+k 

AC-IXl-XRI/RU 
BCMY1-YR  >/R0 
CC  * ( 2 1-ZK  )  /  R  0 

C  DETERMINES  *HE  THE  R  RAY  STRIKES  ENCLOSED  iUREACE 
XSS«XP 
YS  S  * Y  P 
Z  S  S  *  7  P 

CALL  IRSCY(  l*XSS*YSS»ZiSI 

xkr»xr 

YRP  * YR 

7RR-ZR 

CALL  T«SCY(I,XkR,YRR*ZKR  ) 

Pl«XkR-XSS 
P2-YRR-YSS 
P  3" 7  k  K-Z  S  S 
X1«PXV(  IS»1  ) 

Y1»PYV(  I  $  » 1  ) 

Zl-PZVI IS*  1  > 

CALL  TRSCYI  I.Xl.Yl.Zl) 

AI«X1-XSS 
B1-Y1-YSS 
Cl »Z  1-ZSS 
I N»N V ( IS) 

DP  U  1 1  *  1 »  I N 
IC-II+l 

if ( ii.ea.iNi ic-i 

X?»FXV(IS»ICI 
Y2-PYVI  I  S  *  I C  » 

Z  2  *  R  Z  V  I 1  S  *  I C  > 

CALL  TPSC Y(  1»X2* YZ* Z2 I 

A2»X?-XSS 

R2-Y2-YSS 

C2-72-ZSS 

0«fl*UU*C2-fl2*Cll-P?*(Al*C2-ClYA?»+P3*(Al*B?-81*A?) 

IF (u.LT.O.JSETURN 

C  RAY  STRIKES  SLRFACE  CN  ENCLOSED  SIDE  Of  SURFACE  F  OG  E 
A 1  »  A2 
B1-B2 
CI-C2 

U  CONTINUE 

C  RAY  STRIKES  ENCLOSED  SURFACE 
I  SK  •  1 
RETURN 
END 

SUBROUTINE  TRSCYI I*XR»YR*ZR) 

C  CONVERTS  CYLINDRICAL  COORDINATES  INTO  RECTANGULAR  CCORDINATES 

COMNLN7CYL/NCY.NSC(IOI*NSP(IO*lO}»XC(iO)»YC(IO)»ZCUO)f  AE  <101*  BE  (I 
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I  IIU.I  I.J 'JRJJUpi 


ii  i 


00)  *CI  1 lu) *KC( 10) 

CN«U«  l  ll.-CHtl**?| 

AN  *  —  t l (I ) *CE  t  i  )  /CN 

bN»-ee <  i‘» *cf <  it  /cn 

AP-ftE  ( I  )*CN-uN*CE (  1  ) 

BP  »— ( A E  ( I  )*CN-AN*CE  (in 
C  P • At ( 1 ) *  BN- AN ♦ B t  (  1  ) 

R0«Afc  <I>*  tXK-XL(  I)  )  *»t  HIM  YR-YC<  II  )*Ct  U  )*<Z«-ZC  U  >) 

XX«XC  <I)*AF  »  I  >*«i 
YY«YC( l » ♦ « t  11  )*RC 

R»SCRT ( 1 XP-XX)**t +1 Y*-YY ) **2  +  l ZK-ZZ ) **2  > 

A  A  » ( X  R—  X  X  )  /  R 
«8»  (  YR-YY  )  /*> 

CC  »  C  ZR-ZZ  I  /R 

A« AN* AA+bK *BB+CN*CC 

IFU.KT.I.IAM. 

IF ( A. IT  .-I  .  I A ■— I . 

AN &• ACCS  I  A) 

&■  Ap*AA*F  P*«BtCP*CC 
IF ( C.t T .C  .  1 ANG--4NG 
X«»ANG*RC  (  1  ) 

YR-FG 

ZR-R-RC  1 1  > 

RE  TUF \ 

t  n  r 

5 1'  fc  R  l  G  T  1  K  l  IM  tC  US*»  IS*  *R.YP»  ZO  1 
C  DETERMINE  IF  RAY  STRIKES  CONVEX  SUki-ACE 

C0PPCN/CAN/NT»NC'|FC*NF*NA»N0*NV1  10C>*PXV1KG»8)  *PYV<  H’0»8)*PZV(KC 
Of  B  ) 

CUPfCV/CYL/NCYfASCI 10>fNSP(10fl0)f > C 1 1 0 ) » YC U 0 ) » ZC ( 1 C  )  # AE (  1  0  )  » 6  E  1 1 
00)fCt(10)»RC(10) 

CQMf'C.'N/NORN/AXNl  100  )#  AYN(  1001  »  AIM  IOC) 

COrFrN/LINE/AS»BS»CS# XS# YSfZS# ACfBCfCC 
Ck-AXn{IS)*aS  +  AYN1I5)*8S  +  AZNUS)*CS 

IF (Ck.GE.C . IRETukN 

C  RAY  STRIKES  SURFACE  IN  QUTWARO  DIRECTION 
I»1 

IF<  XS  .EO.F  »VUS#  1 ) .  AMZ.rS.cu.PYVlISf  J  ).  AND  .  Z  S  «  E  0 .  PZ  V  (IS#  I ) )  I  •  i  ♦  1 
S«UXNUSI*<PXV<IS»n-X$l*AYN<mMPYvU.»»!)-Y5>*AZN<lS)<'(PZV(IS*l 
01—25 >)/CK 
XR«AS*S*XS 
YR««S*S+YS 
ZR«CS*S+ZS 
I N »N  V ( 1 S  » 

00  10  I  -1»  IN 

ic » r ♦ a 

IF  ( l.Ei.IM  IC-1 
A  1  •  F  X  V  ( iSf I >-XS 

bi-pyv<  is f  n-Ys 
CI-PZVl  ISf  i  )-ZS 
AZ«PX V< IS* 1C  l-XS 
B2 * P Y V ( I S *  I C  I  —  Y S» 

CZ-PZVI  IS»IC)-2S 
P  1  •  jrsf-jr  S 
P2  »YP-YS 
P 3«7P-ZS 

0«Pl*Hl*C?-t»Z*Cl)-P?MAl*C?-tl*AZ>*P3*UlMZ-t'l*A2) 

IF10.LT.C .)RETU-K 
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C  PAY  STRIKES  SURFACE  CN  ENCLOSED  SIDE  OF  SURFACE  EDGE 
10  CONTINUE 

C  RAY  STRIKES  ENCLG'EU  SURFACE 
AC  XN  (  I  S  ) 

BC--AYN(1 S) 

CC«-AZN< I S  ) 

ISK-1 

RETURN 

END 

SUBROUTINE  CONP(ANG.RT.TT) 

C  C  Of  PUT  E  S  INCIDENCE  ANGIE.  REFLECTANCE*  AND  T  RAN  SM  IT  T  ANC  E 
C  NATURAL  LIGHT,  ADDITION  OF  POLARIZATION  COMPONENTS  IGNORED 
C0MMC'NZLINE/AS*tJS*CS*X5,YS#ZS»AC*BC*CC 
C  XN,  INDEX  OF  REFRACTION,  TX,  INTERNAL  TRANSNITT  ANCE 
CATA  XN, TXX/1 . 5, .05/ 

TT-C. 

A«AS*AC+8S*bC+CS*CC 
IF  <  A.LT .C. » A «  —  A 
IMAES(A)  .GT.l.  )  A  •  1  • 

ANG-ACOSU) 

ANGP-  ASIMSIMANGI/XN) 

ca-ccs  <  anc  » 

SA-SINE ANG) 

SI -SCRT ( XN*»?-S A**?  ) 

R0»(((CA-Sl)/(CA  +  Sin+*2+((CAMXN  +  *2t-:i)/(CA*(  XN**2  1  +  S  1 1  )**?)/ 2 
T0«(1.-R0I+CA/CGS  (  ANC-P  I 
CA-CC S(ANGP) 

SA-S  INI  ANC-P  ) 

TX-EXPl— TXX/CA) 

Sl-SCPT (XN**2-SA**2  ) 

PI- << (CA-S1 )/ (CA  +  S1 J >**2  +  < <CA*< XN**2>-fl >✓( CAM XN**2)  +  $1>)  **2)/? 
TI-I1.-RI )*CA/COS(ANG> 

TT  »Tr*TI*Tx/(l.-(»I*TX)**2l 

RT-PC  +  R I*  T  T*T  X 

RETURN 

END 


c 

COBRA  NCCKU0 

FLAT  PLATE 

CANOPY 

DATA, 

MODEL  209  Art-lS  HELICOPTER 

c 

NO. 

VERTICES, 

AND  SliP  E  AC 

f! S  —  BODY,  CONSTRAINT,  AND  TRANSPARENT 

190 

r 

83 

90 

C 

<10 

c 

NO. 

VERTICES 

PEP  SURFACE 

9 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

6 

6 

A 

A 

A 

A 

6 

6 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

5 

5 

5 

5 

c 

SURFACE 

VERTICES  IN 

ROTATIONAL 

nROER 

-  ALL  1ST  VERTICES,  ALL  2ND, 

A 

A 

l 

5 

3 

12 

9 

9 

13 

15 

17 

16 

15 

22 

19 

22 

2  A 

19 

2  5 

20 

29 

38 

35 

A2 

39 

A8 

A3 

A  9 

52 

53 

56 

60 

57 

6A 

61 

65 

7? 

72 

71 

76 

73 

eo 

77 

81 

61 

82 

63 

89 

i 
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C  X-P0S1TILN  All  VtSTICfci 


-6. 

-6. 

6. 

6. 

-17. 

17. 

13. 

-13. 

-2. 

2. 

2. 

-2. 

2. 

2. 

-2. 

-6.25 

-15.5 

-10.907 

-6.25 

10. 093 

5.312 

-16. 

16. 

-10.5 

-10.5 

-10.5 

-10.5 

15.  5 

14. 938 

15.063 

15.063 

-15.125 

15.125 

15.425 

15.675 

13.5 

15.032 

15.032 

13.5 

13. 

11.375 

11.375 

13. 

-11.157 

-12. 782 

-12.782 

-11.157 

-13. 

13. 

13. 

13. 

10. 125 

11.75 

13.344 

11.719 

-12. 

-10.375 

-11.282 

-13.282 

10.782 

12. 282 

17.062 

15.562 

-13. 282 

-11.782 

-15.562 

-16.782 

12.282 

12.282 

-12.282 

-1 2.282 

1C. 657 

12.657 

14.375 

12.375 

-13.25 

-11.232 

-12.375 

-14.375 

10. 125 

-10. 375 

-10.375 

10.094 

12.375 

-11. 313 

13.5 

15.032 

-15.063 

-16.782 

-13.25 

-12.657 

14, 375 

-13.282 

-  1  (  .  6  8  8 

-16.688 

Y-PO^ITT 

c n  ali  vmices 

22.788 

22. 788 

2c . 738 

22.788 

57.288 

57.288 

55.280 

5  5.2  86 

53 


\ 


-12. 

-IP. 

12. 

12. 

-15. 

15. 

2. 

2. 

-2. 

-2. 

-2. 

2. 

-15.5 

15.5 

15.5 

6.25 

-4.969 

-9.969 

6.25 

10.875 

15.313 

5.375 

-5.812 

-16.093 

10.5 

1C.  5 

1C. 5 

10.  5 

-15.5 

-14 . 906 

-14.e75 

-14.675 

15.75 

-15.75 

-15.875 

-15.4 

-13.407 

-15.032 

-15.03? 

-13.407 

11.75 

10.125 

10.125 

11.75 

-10.375 

-12. 

-12. 

-10. 375 

-13. 

-13. 

11.75 

1C.  125 

11.719 

13.344 

-10.375 

-12. 

-13.282 

-11.28? 

12.282 

10. 732 

17.25 

15.75 

-] 1.75 

-13.25 

-16.688 

-17.25 

12.282 

12.2o2 

-12.262 

-12.262 

12.657 

10.657 

12.375 

14.375 

-11.28? 

-13.282 

-14.375 

-12.375 

-11.157 

11.375 

10.657 

-11.262 

-12.375 

11.719 

1 7 • 0b2 

13.625 

12.157 

-13.407 

13.75 

17.25 

16.262 

16.063 

-15.907 

-14.375 

• 

• 

55.288 

55.288 

55.288 

55.288 

5  7 . 2  oe 

57.28b 

57.786 

57.768 

_  * 


57. 78f 

57.786 

67  «  2  8 h 

67,288 

67.28b 

67.28b 

73.188 

73.188 

(>  7  .  2  H  8 

67.288 

73.188 

73.188 

92.768 

42. 766 

95.286 

1C1 . 768 

101. 788 

4tu  288 

98.505 

101.788 

101.853 

99.853 

101.766 

98.249 

90.78 

iul. 875 

110.788 

110.788 

110.25 

11C. 062 

110.26 

110.71 

105. 788 

110.063 

105.760 

107.768 

105. 7( fc 

110.663 

104. 7b6 

107.766 

133.058 

155.125 

151.655 

135.595 

133.038 

155.030 

151.5 

130.566 

50.158 

00.158 

102.908 

155.658 

151.508 

151.508 

105.656 

152.908 

57. 56? 

59.975 

61.575 

69.06? 

57.353 

59.475 

61.575 

56.653 

57. 18  7 

57.187 

57.167 

57.  187 

72.218 

72.218 

72.21b 

72.216 

56.537 

56.637 

56.537 

5b. 537 

72.21b 

72.21b 

72.218 

72.218 

71.375 

71.568 

7?  .2  lb 

73.216 

73.216 

72.218 

72.218 

72.216 

72.718 

72. i  1  3 

110.216 

110.216 

110.218 

110.218 

72.218 

72.216 

72.21  8 

72.216 

110.218 

110.218 

110.218 

110.218 

108. C3 

108.03 

110.568 

110,56“ 

96.537 

98.537 

101.155 

101.  15  5 

106.312 

106.312 

100. 186 

1  0  4 . 1  5  o 

98.537 

48.218 

101. 353 

101.155 

108.28 

108.26 

110.218 

110.216 

106.28 

108.28 

110.218 

108.28 

108.28 

10  8,26 

108. 78 

108. 28 

163.28 

153.28 

153. 2fi 

153.28 

108.312 

108.312 

108.312 

108.312 

IK. 655 

110.654 

11C. 655 

110.605 

51.937 

07.187 

71.568 

71.375 

73.216 

73.216 

106.28 

lOo. 312 

110.655 

110.216 

153.28 

1  5  3 . 1  2  6 

59.062 

61.375 

98.537 

106.03 

73.216 

08.853 

61.312 

98. 218 

10*  .312 

73.26 

108.937 

101.155 

l  55.28 

150.218 

152.28 

100.  16  6 

101.353 

155.376 

106.875 

153.25 

• 

• 

7-1  OS  1  T  It.N  Ul  VERTICES 

66.556 

56.556 

55.556 

0  6  •  5  5  t 

71.246 

55.556 

55.556 

71 .2bo 

61.946 

6 1 . 956 

r>  .5  56 

78.556 

55.556 

55.556 

bb.456 

65.996 

65.95 t 

feb .656 

7  ti  •  5  ■,  6 

75.556 

75.556 

70.556 

78.556 

78.55b 

71.558 

71.556 

71.556 

71.556 

69.556 

69.556 

70.556 

b  5 . 94  6 

85. 958 

75.556 

79.718 

79.956 

83.  353 

83. 730 

79.956 

79.043 

83.687 

►3.358 

7-* .  5  7o 

74.576 

79.53b 

V 3 . 662 

86.537 

78.537 

77.576 

6  0  •  P  7  6 

►3.576 

r>3.6  7fc 

77.576 

■30.076 

63.576 

83.576 

83.076 

37.6b? 

57.187 

35.537 

8  3  .  C  76 

67.625 

87.187 

05.568 

70.821 

70.621 

84.321 

95.696 

105.956 

105.956 

95.696 

89. 321 

77.156 

60.275 

66.2  70 

77.106 

77.25 

6  6. 

66. 

77.25 

76.  25 

7c.  26 

75.25 

75.25 

92.167 

4?.lb7 

40. Ib7 

40.  It7 

76.003 

76.093 

75.093 

75.093 

92.562 

92. 562 

90.562 

90.062 

01.312 

91.626 

42.187 

92.812 

92.562 

42.187 

92.187 

92.187 

90.  187 

90.187 

105.968 

105.968 

103.968 

103.968 

92. 562 

92. 062 

9C.562 

90.562 

105.458 

105.966 

1 0  3 . 9b  3 

103. 9t  8 

105.812 

105.812 

105.812 

105.612 

75.237 

75.237 

76.067 

76.087 

1C5.032 

105.937 

105.812 

105.612 

70.237 

75.087 

75.637 

76.087 

103.812 

105.812 

105.812 

103. 812 

103.812 

105.812 

105.812 

105.812 

105.812 

105.812 

103.81 2 

103. 612 

1C  0.4  06 

105. 9Co 

103.906 

103.906 

105.718 

105.716 

103.718 

103.716 

106. 

106. 

105. 

105. 

76.043 

76.25 

01.625 

9l . 31? 

92 . 662 

92.612 

ICO.fcl? 

105.71" 

106. 

105.96r 

105  .oct 

105.718 

75 .166 

66.370 

75.537 

102.937 

9  C  •  l  b  7 

75.2  5 

66.062 

75.187 

102.937 

90.537 

103.906 

76.137 

88.353 

9  3.75 

1C5.062 

l 33. 812 

70.437 

88.81? 

93.906 

105.124 

• 

• 

CYIING*ICU  DATA 
r\ 

■ 

V 

P  ILOT-t  YE 

60S  I  T I CN 

0. 

138.968 

95.506 

POCKUP  PCSltiro,  ANP  LRIEfTATKlS  DN 

TEST  fLOuR 

0. 

25.75 

-30 

.086 

0. 

-1 . 

0. 

0. 

28 

C. 

1  . 

58. 

-U0.55  36. 

60. 

- b8 . 55 

36. 

36. 

-76.55 

12. 

54 


I 


55 


i 


